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Characterization  of  Stripline  Crossing  by  Transverse  Resonance  Method 

Abstract 

A  method  of  analysis  is  described  for  characterizing  the  discontinuities  made  of 
orthogonally  crossed  two  striplines  on  a  suspended  structure.  The  method  of  analysis  is 
based  on  generalized  transverse  resonant  technique  extended  here  to  4-port  configurations. 
The  technique  is  used  for  determination  of  resonant  structure  at  a  given  frequency  and 
subsequently  the  equivalent  circuit  parameters  of  the  discontinuities. 
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1.  INTRODUCTION  OF  THE  PROBLEM 

Stripline  crossings  of  die  multi-layer  printed  circuitboard  are  commonly  used  in 
digital  circuit  design.  As  the  signal  frequency  gets  higher  due  to  high  speed  processors, 
an  accurate  wave  analysis  of  the  characteristics  of  the  crossing  becomes  important.  In 
addition,  the  crossing  of  strips  on  both  sides  of  the  suspended  substrate  often  appears  in 
microwave  and  millimeter  wave  integrated  circuits  [1].  To  date  little  has  been  reported 
on  the  exact  analysis  of  such  structures. 

The  problem  presented  here  is  to  characterize  the  discontinuities  of  orthogonally 
crossed  two  striplines.  The  structure  to  be  analyzed  is  shown  schematically  in  Fig.  1-1 
along  with  the  coordinate  system.  Two  striplines  are  crossed  orthogonally  on  opposite 
sides  of  the  substrate.  Auxiliary  conducting  planes  are  added  to  convert  the  structure  to  a 
closed  one.  Depending  on  which  ground  plane  (at  z  =  -hi  or  z  =  c)  is  associated  with 
each  strip,  the  structure  may  be  regarded  as  consisting  of  either  inverted  or  suspended 
striplines. 


Region  (2) 
Region  (3) 

Region  (1) 
y=b 


SI  :  Stripline  1 
S2 :  Stripline  2 


x=a 


Fig.  1-1  Structure  for  the  problem 


It  is  assumed  that  each  pair  of  opposing  side  walls  does  not  influence 
electromagnetic  fields  guided  by  the  strip  parallel  to  the  walls  but  only  the  field  guided 
along  the  orthogonal  directions.  This  assumption  is  valid  as  long  as  the  field  remains 
confined  in  the  proximity  of  the  two  striplines.  Thus,  surface  wave  and  radiation 
phenomena  are  excluded.  For  the  sake  of  generality,  two  different  dielectric  layers  with 
dielectric  constant  el  and  e2  are  assumed  in  the  region  (1)  and  (2),  respectively.  The 
auxiliary  walls  are  used  for  field  analysis  purpose.  They  permit  the  structure  to  be 
analyzed  as  a  rectangular  waveguide  discontinuity  problem. 

The  method  for  analysis  is  based  on  a  generalized  "transverse  resonance  technique" 
introduced  for  finline  step  discontinuity  problem  [2].  The  technique  is  extended  here  to  a 
4-port  configuration  treated  in  this  report.  The  method  consists  of  two  parts.  First,  the 

resonant  structure  created  by  auxiliary  walls  is  described  in  terms  of  network 
representation  containing  a  reactive  4-port.  For  a  fixed  resonant  frequency,  we  try  to 
find  as  many  resonator  sizes  as  required  for  extraction  of  4-port  matrix  elements.  The 
second  part  of  analysis  is  a  full- wave  electromagmetic  field  analysis  in  which  the  resonant 
frequency  is  found  as  an  eigenvalue  problem.  A  unique  feature  of  this  second  part  of  the 
analysis  is  to  view  the  problem  as  the  waveguide  scattering  for  the  waves  traveling  in  the 
direction  normal  to  the  subatrate  surface. 


vs 


2.  CIRCUIT  REPRESENTATION 


2*1.  Resonance  Method 


Fig.  2-1  Microstrip  line  with  some  discontinuities 


In  this  section,  the  resonance  method  applied  to  2-port  network  [3]  is  briefly 
described.  Suppose  we  have  a  transmission  line  with  some  discontinuities  as  shown  in 
Fig.  2- 1 .  With  appropriate  reference  planes  for  the  input  -  output  ports,  the 
discontinuities  can  be  represented  as  a  2-port  network  matrix  form  : 


VI 

=  [Z] 

11 

[Z]  = 

Zll 

Z12  ' 

.V2, 

.12. 

* 

_Z21 

Z22  _ 

where  Z1 1,  Z12,  Z21,  Z22  are  unkown  parameters.  In  the  assumption  of  lossless  and 
reciprocal  network,  we  have  only  three  real  unknowns  in  the  matrix.  If  the  two  ports  are 
reactively  terminated,  the  entire  circuit  becomes  a  resonator.  This  is  shown  in  Fig.  2-2 
where  Zoi  (i  =  1,2 )  is  characteristic  impedance,  Bi  is  the  propagation  constant  and  li  is  the 
length  of  the  i  th  transmission  line.  In  Fig.  2-2,  the  resonance  condition  is  expressed  as 
(Zll  +  Z1  )(Z22  +  Z2)-Z122  =  0  .  (2-2) 


reference  plane 


Fig.  2-2  Entire  circuit  as  a  resonator  for  a  2-port  network 


If  the  same  resonant  frequency  is  obtained  with  three  different  pairs  of  Z1  and  Z2,  Eq. 
(2-2)  yields  three  different  equations  and  the  Z  parameters  can  be  solved.  The  reference 
plane  can  be  placed  at  any  position  as  long  as  it  is  on  the  continuous  part  of  the 
transmission  line  and  the  field  disturbance  due  to  the  discontinuities  remains  the  same. 

The  equivalent  circuit  of  a  2-port  network  can  be  commonly  expressed  as  shown  in 
Fig.  2-3. 


II  Z11-Z12  Z22-Z12  12 


Fig.  2-3  Equivalent  circuit  of  a  2-port  network 

Since  the  value  of  Z12  is  solved  with  its  square  value  in  Eq.  (2-2),  the  sign  of  the 
value  must  be  chosen  so  that  the  equivalent  circuit  has  physically  proper  frequency 
performance  in  the  regarding  frequency  region.  If  the  reference  planes  are  placed  close 
to  the  discontinuities,  where  the  fields  are  much  influenced  by  them,  negative  capacitors  or 
inductors  may  appear  in  an  equivalent  circuit  representation. 


2-2.  4-Port  Network 

In  this  section,  a  procedure  for  a  2-port  resonance  method  is  extended  to  a  4-port 
configuration.  The  crossing  between  the  two  suspended  striplines  can  be  represented  as 
a  4-port  network  at  some  reference  planes  sufficiently  far  from  the  discontinuity  region. 
Each  port  is  terminated  with  a  reactance  corresponding  to  the  line  section  between  the 
reference  plane  and  the  auxiliary  wall  as  shown  in  Fig.  2-4. 


Fig.  2-4  4-port  network  for  the  problem 


The  network  equations  for  the  entire  circuit  is  expressed  in  matrix  form  as 


[ V ]  =  £  Z]  [ I ] 

(2-3a) 

[  V  ]  =  -  diag  [  Zi  ]  [  I  ] 

(2- 3b) 

thus 

[  [  Z  ]  +  diag[  Zi  J  ]  [  I  ]  =  0 

(2- 3c) 

where  Zi  ( i  =  1, 2,3,4 )  are  the  terminal  impedance  :  Zi=  j  Zo i  tan  Mi .  The  entire 
structure  is  assumed  to  be  lossless  so  that  the  Z  parameters  of  the  4-port  network  are 
imaginary  and  the  resonant  frequency  is  real.  The  resonant  frequency  is  obtained  from 
the  condition  that  the  voltages  and  currents  are  non-trivial  in  the  absence  of  sources.  From 
Eq.  (2-3c),  this  condition  is 

det  II  [  Z  ]  +  diag  [  Zi  ]  II  =  0  (2-4) 

The  values  of  Zi' s  can  be  specified  once  the  distance  to  the  wall  is  fixed  while  the  Z 
parameters  are  to  be  determined.  The  impedance  matrix  of  a  reciprocal  4-port  lossless 


6 

network  possesses  in  general  10  independent  imaginary  parameters.  In  the  present  case, 
however,  because  of  the  symmetry  of  the  structure,  only  five  parameters  are  needed  to 
characterize  the  Z  matrix : 

VI  1  f  Zll  Z12  Z13  Z13 1  f  II  ' 

V2  Z12Z11 Z13  Z13  12 

(2-5) 

V3  Z13  Z13  Z33  Z34  13 

V4  J  [zi3  Z13  Z34  Z33  J  l  14  „ 

From  Eq.  (2-5),  the  resonant  condition  of  Eq.  (2-4)  can  be  written  in  the  following  form  : 

{  (Z11+Z1)(Z11+Z2)-Z122  }  {  (Z33+Z3)(Z33+Z4)-Z342  } 

-  4  Z132{  Z1  l-Z12+(Zl+Z2)/2  }  {  Z33-Z34+(Z3+Z4)/2  }  =  0  (2-6) 

The  derivation  of  the  equation  is  detailed  in  Appendix  A.  It  is  observed  that,  in  the  limit 
as  Z13  equals  zero,  Eq.  (2-6)  reduces  to 

(Z11+Z1)(Z11+Z2)-Z122  =0  or  (2-7a) 

( Z33  +  Z3 )  ( Z33  +  Z4 )  -  Z342  =  0  .  (2-7b) 

As  expected,  when  the  two  striplines  are  uncoupled,  the  resonant  condition  splits 
into  those  of  the  individual  lines.  The  unknown  Z  parameters  as  well  as  terminal 
impedances  Zi’s  are  frequency  dependent.  Eq.  (2-6)  can  be  regarded  as  a  function  of  to 
equated  to  zero.  The  Z  parameters  in  Eq.  (2-6)  may  be  solved  if  the  five  different 
equations  are  derived  corresponding  to  five  different  set  of  Zfs  at  t  ie  same  resonant 
frequency. 

Now  we  can  show  that  by  properly  choosing  the  terminal  impedance  Zi’s,  the 
resonant  conditions  are  simplified  so  that  the  problem  is  solved  analytically.  First,  let  us 

choose  the  terminal  impedances  in  a  symmetrical  way,  i.e.,  11=12  and  13=14  so  that 

Z\  =72  (2- 8a) 

Z3  =  Z4  .  (2-8b) 

The  above  conditions  correspond  to  the  side  walls  located  symmetrically  with  respect  to  the 
discontinuities.  If  Eqs.  (2-8)  are  applied,  Eq.  (2-6)  can  be  factorized  in  the  form 


(Zl 1+Z1-Z12)  (Z33+Z3-Z34)  {(Zl  l+Z12+Zl)(Z33+Z34+Z3)-4  Z13^}  =  0  .  (2-9a 

Thus  Zll  +  Zl  -Z12  =0  or  (2-9b 

Z33  +  Z3-Z34  =0  or  (2-9 

(Zll  +  Z12  +  Z1  )(Z33  +  Z34  +  Z3)-4Z132  =  0  .  (2-9 

With  each  factor  in  Eq.  (2-9a)  equated  to  zero,  the  eigenvalues  for  the  matrix  in  Eq.  (2- 3c) 
are  obtained ;  each  eigenvalue  o>  is  then  the  resonant  frequency  under  the  condition  of  the 
corresponding  eigenvector  of  the  currents.  When  the  first  factor  in  Eq.  (2-9a)  is  equated 
to  zero,  i.e.  Eq.(2-9b)  is  satisfied,  the  eigenvector  of  Eq.  (2-4)  is  easily  found  as  detailed 


in  Appendix  B : 


I  1  =  -12 
13=14  =  0  . 


(2- 10a) 
(2- 10b) 


This  condition  corresponds  to  an  odd  resonance  of  the  structure  shown  in  Fig.  2-5(a). 

The  structure  behaves  as  if  an  electric  wall  is  placed  symmetrically  along  at  the  center  of 
stripline  2.  For  the  given  resonant  frequency,  the  required  resonance  condition  provides 
the  quantity  Zl  1-Z12  for  a  given  value  of  Zl.  13  =  14  =  0  ( therefore  V3  =  V4  =  0  unless 
Zi  is  infinite )  does  not  imply  that  the  electromagnetic  field  is  zero  along  the  strip  2,  but 
only  that  the  electromagnetic  field  of  the  dominant  mode  is  zero ;  it  is  anticipated  that  the 
field  is  localized  at  the  crossing  region  along  the  strip  2  as  illustrated  in  Fig.  2-6(b). 


n  Z11-Z12  Z11-Z12  u 


13  Z33-Z34  Z33-Z34  I3 


Zl  1  -Zl  2 


(a)  stripline  1 


Z33-Z34 


(b)  stripline  2 


Fig.  2-5  Equivalent  circuits  for  an  odd  resonance 


(a)  Cross  sectional  view  (b)  Side  view 

Fig.  2-6  Longitudinal  current  flow  on  stripline  2  with  E-wall 


Similary,  when  Eq.  (2-9c)  is  satisfied,  an  odd  resonance  of  stripline  2  is  obtained 

11=  12  =  0  (2- 11  a) 

13  =  -14  .  (2- lib) 

The  resonant  circuit  is  expressed  in  Fig.  2-5(b).  It  is  interesting  to  note  that  the 
individual  parameters  of  the  pairs  of  ports  1-2  and  3-4  are  uncoupled  in  the  odd  resonance 

modes.  Finally,  from  Eq.  (2-9d),  the  eigenvector  for  an  even  resonance  is  obtained  : 

11=  12  (2- 12a) 

13  =  14  (2- 12b) 

Substitution  of  this  condition  into  Eq.  (2-5)  yields  the  2-port  network  matrix  equation 


VI 

-V3 


Z11+Z12  2  Z13  " 

• 

11 

.  2Z13  Z33+Z34^ 

.13  . 

(2-13) 


whose  equivalent  circuit  is  illustrated  in  Fig.  2-7  together  with  terminal  impedances. 


Z11'  =  Z11  +Z12 
Z12'  =  2  Z13 
Z22'  =  Z33  +  Z34 


Fig.  2-7  Equivalent  circuit  for  an  even  resonance 
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From  the  condition  of  Eqs.  (2-12),  it  is  postulated  that  the  even  resonance  structure  has 


two  magnetic  walls  at  each  center  of  the  striplines.  The  use  of  symmetry,  therefore,  has 


reduced  the  4-port  network  problem  to  that  of  a  2-port.  For  a  given  resonant  frequency, 


three  different  pairs  of  Z1  and  Z3  ( with  2,1-72  and  Z3  =  Z4 )  are  used  to  compute  three 


quantities  Z1 1  ',Z12'  and  Z22'  which  denote  the  elements  of  the  matrix  in  Eq.  (2-13) : 

Zll'  =  Zll  +  Z12  (2- 14a) 


Z22'  =  Z33  +  Z34 


Z12'  =  2  Z13 


(2- 14b) 
(2- 14c) 


Combining  the  results  with  those  for  the  two  structures  in  Fig.  2-5,  we  obtain  all  five  Z 


parameters.  For  the  procedure  illustrated  above,  we  must  know  the  propagation 


constants  of  the  two  isolated  striplines,  i.  e.,  B1  and  B3.  The  quantities  are  necessary  to 


obtain  Z1  and  Z3  : 


zl  =  j  tan  Bl(  a  -w2  )/2 
z3  =  j  tan  B3(  b  -wl  )/2 


(2- 15a) 
(2- 15b) 


where  the  impedances  are  normalized.  These  expressions  are  obtained  for  a  specific 


choice  of  the  reference  planes  as  shown  in  Fig.  2-8.  The  field  analysis  can  be  also  used 


to  determine  B1  and  B3.  This  can  be  done  simply  by  setting  B1  =  jc/a  and  B3  =  rc/b  at  a 


specified  resonant  frequency  for  the  isolated  stripline. 


reference  planes 


wl  si 
B1 


.  2 


a-w2 
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Fig.  2-8  A  choice  of  reference  planes  (top  view) 
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2-3.  Equivalent  Circuit 

It  is  important  to  find  an  equivalent  circuit  representation  of  the  4-port  structure  so 
that  it  is  most  convenient  for  the  present  analysis.  Intuitive  insights  are  needed  to  figure 
out  a  good  expression  which  well  represents  a  structure  in  a  broad  frequency  range.  A 
basic  idea  for  choosing  proper  elements  is  to  use  the  knowlege  of  field  distributions.  If 
the  E  field  originating  from  a  stripline  increases  in  the  proximity  of  discontinuities,  it  may 
be  represented  by  a  parallel  capacitive  element  When  the  H  field  circling  around  a 
stripline  changes,  it  may  be  represented  by  a  series  inductive  element. 


Fig.  2-9  Equivalent  circuit  for  the  problem 


For  the  structure  in  Fig.  1-1,  one  possible  equivalent  circuit  representation  could  be 
chosen  as  shown  in  Fig.  2-9.  The  circuit  already  takes  the  symmetry  properties  into 
account  Zc  is  used  to  represent  the  coupling  capacitance  between  the  two  strips. 

Zp+Zp  and  Zq+Zq  represent  the  inductances  associated  with  the  two  stripline  sections, 
while  Za  and  Zb  represent  the  two  strip-to-ground  capacitances.  The  Z  matrix  for  the 


circuit  in  Fig.  2-9  is  given  as 


Zll  Z12Z13Z13 


Z12Z11  Z13Z13 
Z13Z13Z33  Z34 


(2- 1 6a) 


where 


Z13Z13Z34  Z33 


Zll  =Zp  +  Za(Zb  +  Zc)/(Za  +  Zb  +  Zc) 
Z12  =  Za  ( Zb  +  Zc )/( Za  +  Zb  +  Zc ) 

Z13  =  Za  Zb/(  Za  +  Zb  +  Zc ) 

Z33  =  Zq  +  Zb(Zc  +  Za)/(Za  +  Zb  +  Zc) 
Z34  =  Zb  ( Zc  +  Za )/( Za  +  Zb  +  Zc ) 

From  Eqs.  (2-16),  each  element  in  the  equivalent  circuit  is  obtained : 

Za  =  ( Z132  -  Z12  Z34 )/( Z13  -  Z34 ) 

Zb  =  ( Z132  -  Z12  Z34  )/( Z13  -  Z12 ) 

Zc  =  -  ( Z132  -  Z12  Z34  )/Z13 
Zp  =  Zll  -Z12 
Zq  =  Z33  -  Z34 


(2- 16b) 
(2- 16c) 
(2-16d) 
(2-16e) 
(2-160 

(2- 17a) 
(2- 17b) 
(2- 17c) 
(2-17d) 
(2-17e) 


The  derivation  of  Eqs.  (2- 16)  and  (2-17)  are  detailed  in  Appendix  C.  As  computation  for 
the  Z  parameters,  accordingly  the  parameters  in  Eqs.  (2-17),  it  is  convenient  to  use 
normalized  impedance  parameters  with  respect  to  the  characteristic  impedance  of  each 
stripline.  Thus,  the  matrix  in  Eq.  (2-3)  may  be  modified  as  follows : 

fvl  I  fzll  z!2  zl3  zl3l  f  il  1 


zl2  zll  zl3  zl3 
zl3  zl3  z33  z34 
zl3  zl3  z34  z33 


(2- 18a) 


where 


[  v  ] :  normalized  voltage  vector 
[  i  ] :  normalized  current  vector 

[  v  ]  =  diag  [  VZoTr1!  V  ] 
[  i  ]  =  diag  [  Vz3T ]  [  I  ] 


zmn  = 


/Zom  VZon 


m  =  1, 2,  3, 4 
n  =  1,2,  3, 4 


(2- 18b) 
(2- 18c) 

(2- 1 8d) 


Zo i :  Stripline  characteristic  impedance  at  the  i-th  port 
The  normalized  terminal  impedance  z/'s  are  given  by  Eqs.  (2-15).  Small  letters  z,  v,  i 
for  Z,  V,  I  denote  normalized  quantities.  Accordingly,  the  equations  in  Eqs.  (2-8),  (2-9) 

and  (2-13)  are  replaced  by  normalized  ones,  -  replacing  Z's  by  z's  - ,  and  solved  with 


respect  to  normalized  parameters.  Finally,  equivalent  circuit  parameters  in  Eq.  (2-17)  are 
solved  using  the  relation  of  Eq.  (2-18d). 

The  S  matrix  for  a  4-port  network  is  obtained  form  the  relation 

[S]  =  [[z]  +  [U]]'1[[z]-[U]]  (2-19) 

where  [  z  ]  is  a  normalized  impedance  matrix  and  [  U  ]  is  a  unit  matrix.  The  matrix 
relations  are  detailed  in  Appendix  D. 
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3.  FIELD  ANALYSIS 

3-1.  Modes  in  a  Rectangular  Waveguide 


y 

b 


0 


cross  sectional  area 


E  or  M  wall 


a  x 


Fig.3-1  Rectangular  waveguide 


Fig.  3-1  shows  a  rectangular  waveguide  structure  with  the  reference  coordinate. 
With  the  scalar  potentials,  the  general  field  expressions  for  TE  to  z  and  TM  to  z  modes  are 
as  follows  [4] : 


TE 

Et  =  -VtVxz 


H  t  =  -  Vt  \jr 
zA 

kc2 

Hz- - y 


TM 

H  t  =  Vt  <P  x  z 


where 


E  t  = 

_j£_ 

yA 

kc2  ^ 

E  z  = 

yA  ' 

Vtq> 


<P 


(3-1) 


kc2  =  ( -Q22L  )2  +  (  JlZL  )2 
a  b 

m,  n :  see  Eqs.  (3-2) 

B2  =  k2  -  kc2 

yA=jcoe,  zA=jto^i,  k2  =  -yAzA  =  co2en 
e :  permitivity  p. :  permeability 
E  t  =  x  Ex  +  y  Ey 
H  t  =  x  Hx  +  y  Hy 

Vt  =  x-^-+y-^-  x,  y,  z :  unit  vector  in  x,  y,  z  direction 
ox  oy 


where  z  dependanee  e‘J“z  is  assumed  and  subscript  t  denotes  the  transverse  components  of 
a  quantity.  Assuming  an  electric  wall  or  a  magnetic  wall  at  each  center  of  the  wave 
guide,  one  can  express  scalar  potentials  including  conductive  guide  boundary  conditions  as 


follows : 


Vmn  »  Pmn  cos  ( JESS- )  cos  ( -Ep- ) 
a  b 

'Pmn  =  Pmn  sin  (  EUK. )  sin  (  Ejty  ) 
a  b 

Pmn  =  yl  5m  5n  J_ 
ab  kc 


m,n  =  0,1, 2,3,-- 
m,n  =  1,2,3,  - 


Si  =  {  1  ,  =  0 
1  2  i*0 


(3-2a) 

(3-2b) 


where  subscripts  mn  corresponds  to  the  order  of  a  harmonics,  and  Pmn  is  a  coefficient  to 


normalize  the  norm  of  scalar  potential  vector  functions  so  as 

/  I  Vt  Vmn  £  dxdy  =  1 
S 

J  I  Vt  Vmn  I ^  dxdy  =  1 . 

S 


(3-3a) 


(3-3b) 


Even  symmetry  solutions  with  respect  to  the  x  or  y  axis  are  obtained  for  m  or  n  being  an 
odd  integer.  Similarly,  odd  solutions  are  obtained  for  m  or  n  being  an  even  integer. 


3-2.  Field  Expansions  inside  the  Cavity 


y  -hi 


x  t 


(a)  Region  (1)  -hl£z£0  (b)  Region  (3)  O^z^t  (c)  Region  (2)  t<z^c 


Fig.  3-2  Subregions  for  field  analysis 


The  resonant  structure  of  Fig.  1-1  can  be  subdivided  into  three  homogeneous 
regions  as  shown  in  Fig.  3-2.  In  each  region,  we  can  expand  the  electromagnetic  field  in 
terms  of  TE-to-z  and  TM-to-z  waveguide  modes.  Incorporating  the  boundary  conditions 
for  z=-hl  and  z=c,  the  following  expressions  can  be  assumed  : 
region  (1) :  -hl£z£0 

MN 

Et(1)=  -X  [  Amn(1)Vt  Vmn  x  z  +  Bmn(1)(pmn(1VyAi)  Vt  ^mn  ]  sin  Pmn(1)(z+hl)  (3-4a) 
mn 
MN 

Ht(1)=  X  [  Amn(1)(Pmn(1)/zAj)  Vt  Vmn  +  Bmn(1)  Vt  ^mn  x  z]  cos  Pmn(l)(z+hl)  (3-4b) 
mn 

Pmn(1  >=  ( k  j 2  -  kc2  )  ^  (3-4c) 


region  (2) :  t  <.  z  £  c 

MN 

Et(2>=  -X  [  Amn(2)Vt  Vmn  x  z  +  Bmn^2^(pmn^2VyA2)  Vt  ^mn  ]  sin  Pmn^2)(z-c)  (3-5a) 

mn 
MN 

Ht(2)=  X  [  Amn(2)(Pmn(2)/zA2)  Vt  Vmn  +  Bmn(2)  Vt  'Pran  x  z]  cos  pmn(2)(z-c>  (3-5b) 

mn 

pmn^2^=  (  k2^  -  kc2  (3-5c) 


region  (3) :  0  <,  z  <  t 

MN 

Et(3)=  -X  [  {  Cmn  sin  pmn(3)z  +  Dmn  sin  Pmn(3)(z-t)  }  Vt  Vmn  x  z 
mn 

+  {  Fmn  sin  Pmn(3)z  +  Gmn  sin  pmn(3)(z-t) }  (pmn(2)/yA3)  Vt  ‘Pmn  (3-6a) 

MN 

Ht(3)=  X  [  {Cmn  cos  Pmn(3)z+Dmn  cos  Pmn(3)(z-t) }  (Pmn(3VzA3>  Vt  Vmn 
mn 

+  {  Fmn  cos  pmn<3)z  +  Gmn  cos  pmn^3Yz-t) }  Vt  Vmn  x  z  ]  (3-6b) 


Pmn(3)=  ( k3^  -  kc2  ) ^  (3-6c) 


where  Amn(1),  Amn(2),  Bmn(1),  Bmn(2),  Cmn,  Dmn,  Fmn,  and  Gmn  are  unknown 
coefficients,  and  superscripts  in  parentheses  refer  to  respective  regions.  The  boundary 
conditions  at  each  interface  are 


Interface  (I) :  at  z  =  0 


Etj  =  Et3 

on  S 

(3-7a) 

Et|  =  Et3  =  0 

on  SI 

(3-7b) 

Ht3  -  Ht!  =  {  AHt^ 

on  SI 

0 

on  S-Sl 

(3-7c) 

Interface  (II) :  at  z  =  t 

Et2  —  Et3 

on  S 

(3- 8a) 

Et2  =  Et3  =  0 

on  S2 

(3-8b) 

Ht2  -  Ht3  =  {  AHt(2) 

on  S2 

0 

on  S-S2 

(3-8c) 

where  AHtw  ( i  =  1,2 )  are  unknown  functions  of  x,  y. 

These  functions  are  expanded 

terms  of  a  set  of  known  orthogonal  vector  functions  XyW  defined  over  Si  with  unkown 

coefficients  Pvw  (see  next  section) : 

AHt(i)  =  lPv<»)Xv(1) 

(3-9a) 

V 

AHt(2)  =  Z  Pv(2>  Xy<2) 

(3-9b) 

V 

Applying  the  boundary  conditions  of  Eqs.  (3-7),  (3-8)  to  Eqs.  (3-4),  (3-5),  (3-6) 
and  making  use  of  the  orthogonal  properties  of  Vt^mn  and  Vt^mn,  we  obtain  a 
homogeneous  system  of  equations  in  terms  of  unknown  coefficient  Pv  : 

ViM  VoM 


M 


/im  voiyi 

y  Pv^)  Upv(11)  +  X  Pv(2)Unv(I2)  =  0  p=l,2, 3,  vj 

1  2 

t/iM  VoM 

^  Pv^Upv*21*  +  ^  Pv<2>  Upv(22)  =  0  n=  1,2,3,  v2M 


where 


MN 


c 


(■ 


U  is  A  amn(l^  ^mn,v^  -  (fttyy A.)  bmn(i^  ®mn,pw  ®mn,v^ }  (' 


=  /  XpW  •  VtVmn  ds 
Si 


/=  1,2 


®mn,p^  =  |  XjjW  .  Vt^mn  x  z  ds  i  =  1,2 
Si 


(3-10d) 


(3-10e) 


amn< 1 »  =  Ka  ( 1/S  x )  ( B<2>Ct2  +  B(3)Ct3  )  (3- 100 

amn(12)  =  Ka  ( 1/  (S^) )  B<3>  (3-l0g) 

amn(21)  =  -  Ka  ( 1/  (S2S3) )  B<3>  (3-  lOh) 

amn(22)  =  -  Ka  ( 1/S2 )  ( B<1>Ct1  +  B<3>Ct3  )  (3-10i) 

Ka  =  zA  /  {B(3)2-B(3>Ct3(B<2>Ct2+B<1>Ct1)-B(1>B(2>Ct1Ct2  } 
bmn(11)  =  Kb  ( 1/S1 )  (B(3)Ct2+(e3/e2)B(2)Q3)  (3-10j) 

bmn<l2>  =  Kb  ( 1/  (S1S3) )  (eye2)  B(2)  (3-  10k) 

bmn(20  =  -  Kb  (  1/  (S2S3) )  (e^)  B<0  (3-101) 

bmn^22)  =  -  Kb  ( 1/S2 )  (B<3>Ct1+(e3/e1)B<1)Q3)  (3-10m) 

Kb  =  B<3>/  [  (e3/e1)(e3/e2}B(1)B(2)-B(3)Ct3{(e3/e1)B(1)Ct2+(e3/E2)B(2)Ct1 }  -B<3>2Ct1Ct2  ] 


(3-10g) 

(3-10h) 

(3-10i) 


(3-10j) 
(3- 10k) 
(3-101) 
(3-  10m) 


B(i)  =  Bmn© 


i  =  1, 2,  3 


Cj,  C2,  C3  =  cos  B(1>hl,  cos  B(2)h2,  cos  B(3)t 

Sj,  S2,  S3  =  sin  B(1)hl,  sin  B(2)h2,  sin  B(3)t 

Ctj,  Ct2,  Ct3  =  Cj/Sj,  C2/S2,  C3/S3 

zA  =  j  0)  nQ  (i0  :  free  space  permeability  =  47txlO'2[H/m] 

The  derivations  of  Eqs.  (3- 10)  are  detailed  in  Appendix  E.  The  condition  for  non  trivial 
solution  determines  the  characteristic  equation  of  the  given  structure.  This  equation  may 


be  regarded  as  a  function  of  to,  /1 , 13  equated  to  zero  : 

f  ( to,  /l,/3)  =  0 


(3-11) 


For  given  value  of  0)  =  tor,  Eq.  (3-1 1)  can  be  solved  to  evaluate  the  different  pairs  of  /I 
and  13  giving  rise  to  the  same  resonant  frequency  cor  These  values  of  /I  and  13  can  be 
used  for  computing  the  discontinuity  parameters  discussed  in  the  previous  section. 


3-3.  H  Field  (  Current )  Expansion  on  the  Strips 

The  H  field  on  the  strip  are  expanded  in  terms  of  known  basis  functons  Zy  with 
unknown  coefficients  Pv.  Actually  what  will  be  expanded  is  the  H  field  discontinuity 


which  is  equal  to  Jt  x  z,  viz. 


vmax 

AHt  =  £  Pv  Xy  =  Jt  x  z 
v 


(3-12) 


The  basis  functions  are  chosen  in  such  a  way  that  the  field  is  non-zero  only  on  the  strip. 
Additional  boundary  condition  must  be  satisfied  at  the  strip  ends  where  the  striplines  are 
terminated  with  electric  walls.  The  field  may  be  expanded  in  terms  of  the  basis  functions 
of  two  different  types.  One  of  them  is  harmonic  basis  functions  and  another  is  singular 
basis  functions.  The  singular  behavior  of  the  magnetic  field  component  normal  to  the 
stripline  edges  is  incorporated  in  the  singular  basis  functions.  They  are  therefore 
expected  to  provide  a  faster  numerical  convergence.  For  the  H  field  expansion  on  S 1 
[  see  Fig.  3-3(a)  ],  the  following  set  is  employed  : 


v  fn  .  rrtx  sjc  .  (b-wl) 

*x,rsU)  =  sm -  sin  [ - {y - )] 

a  wl  2 


Zy,rs(1)  =  cos 


COS  [  — {  y  .  }] 

wl _ 2 

V,  ,  (y-b/2)  '2 


(3- 13a) 


(3- 13b) 


(a)  at  interface  (I) 


(b)  at  interface  (II) 


Fig.  3-3  Striplines  at  each  interface 


■\-V  .  /■  .*<  .->■  -v  -*■ .•  .  -  ‘ 


where  subscript  rs  is  used  instead  of  v  for  representing  the  variations  in  the  two  orthogonal 
directions  (see  Appendix  F).  Similarly,  for  the  H  field  on  S2  in  Fig.  3-3(b),  the 
expansions  are 


Xx,rs<2)  = 
Xy,rs(2>  = 


cos  [  QL.  {  x  -  }  ] 

_ )x2 _ 2 _ 


l 


.  (x-a/2)  2 

1  w2/2  ' 


cos 


S7ty 


.  r  nc  (a-w2) 

sm  [ { x - 

w2  2 


sin 


sity 


(3- 14a) 

(3- 14b) 


Depending  on  whether  an  electric  or  magnetic  wall  is  assumed  to  be  at  (A)  and  (B)  in  Fig. 
3-3(a),  the  numbers  r  and  s  in  Eqs.  (3-13)  are  classified  to  be  odd  or  even  as  illustrated  in 
Table  3-1. 


Table  3- 1  rs-table  with  respect  to  E  M  walls 


(a)  at  interface  (I) 


(b)  at  interface  (II) 


4.  COMPUTED  RESULTS 


In  accordance  with  the  technique  described  above,  the  electromagnetic  fields  in  each 
region  are  expressed  in  termes  of  the  series  expansions.  In  the  numerical  computation, 
only  a  finite  number  of  terms  can  be  retained  in  the  series  expansions.  In  order  to  obtain 
a  proper  convergent  behavior  of  the  solution,  the  number  of  the  terms  was  chosen  in  such 
a  way  that  the  highest  spatial  frequencies  of  the  electromagnetic  field  arc  approximately  the 
same  in  the  transverse  directions  (x,y).  The  current  on  a  stripline  is  expanded  in  such  a 
way  that  the  distributions  of  each  component  (Jx,  Jy)  in  the  transverse  direction  are 
expressed  by  only  the  first  term  and  that  the  distributions  in  the  longitudinal  direction  by  a 
set  of  terms,  the  number  of  which  are  the  same  on  each  stripline. 

The  method  is  first  tested  by  computing  the  resonant  frequency  of  the  isolated 
stripline  with  certain  structural  parameters.  Fig.  4-1  shows  the  convergence  of  the 
resonant  frequency  with  respect  to  the  number  of  the  terms  of  the  field  expasions  in  the 
homogeneous  region  of  Fig.  3-2.  The  number  MN  was  chosen  to  be  the  same  in  both  x 
and  y  directions.  The  results  exhibit  very  good  agreement  with  those  by  the  SDA  [5]  in 
which  the  resonant  frequency  was  calculated  from  the  equation 

=  (4-1) 

a  u 

where  c  :  speed  of  the  light 
a :  stripline  length 

B :  propagtion  constant  in  proximity  of  fo 
Bo :  »  in  free  space 

because  computational  the  field  was  represented  by  the  same  number  of  harmonics  in 
both  this  method  and  SDA. 

Fig.  4-2  shows  the  convergence  of  the  resonant  frequency  of  crossing  stripline 
structure  with  respect  to  the  number  v  of  the  current  expansion  terms  in  the  longitudinal 
direction.  At  beginning,  we  expected  that  the  number  of  the  terms  would  be  required  as 


NS 


a  =  30  mm 
b  =  52  mm 


MN  =  100 


Number  of  Basis  Functions  V 

hl=H2  =  5mm,  t  =  1mm 
wl  =  w2  =  1  mm,  er  =  3.8 

Fig.  4-2  Convergence  of  the  resonant  frequency  of  crossing  striplines 

with  respect  to  the  number  of  basis  functions 
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large  as  that  of  the  electromagnetic  field  so  that  the  longitudinal  current  would  have  the 
same  highest  spatial  frequency  and  the  results  would  have  a  good  convergence.  From 
Fig.  4-2,  however,  we  could  observe  a  fast  convergence  characteristic  with  a  relatively 
small  number  of  current  expansions  if  the  electromagnetic  field  is  well  expressed  by  an 
adequate  number  of  harmonics. 

Fig.  4-3  shows  the  parameters  of  the  structure  to  be  computed.  The  substrate 
was  placed  symmetrically  between  the  top  and  bottom  planes  with  both  strips  having  the 
same  widths.  Hence,  the  impedance  matrix  representation  of  the  discontinuities  has 
zl  l=z33  and  zl2=z34,  i.e.  unknown  parameters  number  is  three:  zl  l,zl2,  zl3.  The 
computation  was  done  at  the  three  different  frequencies:  0.5, 1, 2GHz.  The  approximate 
lengths  of  striplines,  namely  the  size  of  cavity  at  each  frequency,  are  also  illustrated  in  Fig. 
4-3.  The  reference  planes  are  defined  in  the  same  mannar  as  illustrated  in  Fig.  2-8. 

Table  4-1  shows  the  element  values  for  the  equivalent  circuit  of  the  discontinuities 
shown  in  Fig.  4-4.  Fig.  4-5  shows  the  corresponding  S  parameters  of  the 
discontinuities.  Note  that  the  elements  for  Za  and  Zb  (=Za)  are  capacitors  with  negative 
value.  This  is  acceptable  because  they  compensate  the  parallel  distributed  capacitance  for 
the  isolated  stripline  with  the  absence  of  the  other  stripline  as  follows.  Assuming  the 
absence  of  stripline  2  by  making  the  port  3  and  4  open,  the  structure  is  nearly  an  isolated 
stripline.  The  section  between  portl  and  2,  therefore,  would  have  an  equivalent  parallel 
capacitor  which  is  coresponding  to  the  distributed  capacitor  of  the  stripline.  When  the 
section  is  very  short,  the  capacitance  of  the  equivalent  capacitor  may  be  negative  so  that  it 
would  compensate  the  fringing  capacitance  due  to  the  open  ends  at  port  3  and  4  whose 
value  could  be  larger  than  that  of  the  distributed  capacitor  as  a  simple  stripline  section. 

Zc  was  also  calculated  with  hl=10mm,  h2=  100mm,  t=5mm,  wl=w2=0.4mm, 
e^lat  2  GHz  simulating  the  crossing  model  of  the  work  by  Giri,  et  al  [6].  The  result  for 
the  value  of  Cc  was  0.082pF,  while  the  value  of  the  coupling  capacitor  Cm  was  estimated 
to  be  0.1 5pF  by  using  the  equivalent  radius. 


a  =  230  -  240 
b  =  200  -  210 


Unit  [mm] 


@  0.5  GHz 

a=  100-  120 

b  =>  90  -  100  @  1  GHz 

a  =  30  -  40 

b  -  50  -  52  @2  GHz 

Fig.  4-3  Structural  parameters  for  computation 


Fig.  4-6(a),(b)  show  the  longitudinal  and  transverse  components  of  the  current 
densities  at  the  center  of  the  strip  at  the  three  different  frequencies.  In  Fig.  4-6(a),  it  is 
observed  that  each  figure  was  of  perturbated  cosine  form  and  had  larger  absolute  amplitude 
than  cosine  in  the  proximity  of  the  center.  This  feature  is  anticipated  from  the  insight  that 
the  stripline  whose  length  is  slightly  shorter  than  a  half  wavelength  turn  out  to  be  a 
resonator  with  capacitive  loading  at  the  center  of  the  stripline.  The  longitudinal  current, 
hence,  may  be  represented  by  the  combination  of  a  cosine  and  an  additional  polynomial 
functions  [7]  so  that  the  computation  time  could  be  reduced.  As  to  the  transverse  current 
shown  in  Fig.  4-6(b),  the  figures  were  not  readily  characterized.  This  implies  that  the 
current  function  would  require  a  number  of  terms  of  basis  functions;  the  current  function 
could  not  be  easily  replaced  by  few  terms  of  special  functions. 

Fig.  4-7(a)(b)(c)  show  the  computed  results  for  different  values  of  the  stripline 
width  and  the  substrate  thickness. 

5.  CONCLUSIONS 

A  method  of  analysis  has  been  described  for  characterizing  the  discontinuities  of 
two  crossed  striplines.  The  method  is  based  on  a  generalized  transverse  resonance 
technique  for  computing  the  resonant  frequency  of  a  resonator  created  by  enclosing  the 
crossing  with  auxiliary  perfectly  conducting  walls.  This  resonator  problem  is  analyzed  as 
the  waveguide  scattering  for  waves  traveling  in  the  direction  normal  to  the  substrate 
surface.  For  a  specified  frequency,  resonant  structures  are  found  by  adjusting  the 
lengths  of  the  strips  and  hence  the  resonator  size.  These  structures  are  used  for  deriving 
the  equivalent  circuit  parameters  characterizing  the  discontinuity. 

This  method  can  also  be  applied  for  the  characterization  of  stripline-slotline 


transition. 


APPENDIX 


A.  DERIVATION  OF  EQ.  (2-6) 
det  II  [  Z  ]  +  diag  [  Zi  ]  II 

=  fzil+Zl  Z12  Z13  Z13 
Z12  Z11+Z2  Z13  Z13 
Z13  Z13  Z33+Z3  Z34 
Z13  Z13  Z34  Z33+Z4 


=  Z11+Z1-Z12  Z12-Z11-Z2  0 


Subtract  row  2  from  row  1 
Subtract  row  3  from  row  4 


Z12  Z11+Z2 


Z13  Z13 

0  0 


Z33+Z3  Z34 


Subtract  column  2 
from  column  1 

Subtract  column  3 
from  column  4 


Z34-Z33-Z3  Z33+Z4-Z34 


=  2Z11-2Z12+Z1+Z2  Z12-Z11-Z2 


Z12-Z11-Z2 


Z11+Z2 


Z33+Z3 


Z34-Z33-Z3 


L  0  0  Z34-Z33-Z3  2Z33-2Z34+Z3+Z4 

=  {  (  2Z1 1-2Z12+Z1+Z2  )  (  Zl  1+Z2  )  -  ( Z12-Z1 1-Z2  )2  ) 

x  (  ( 2Z33-2Z34+Z3+Z4 )  (  Z33+Z3  )  -  ( Z34-Z33-Z3  )2  } 

-  4  Z132  {  Zll-Z12+(Zl+Z2)/2  }  {  Z33-Z34+(Z3+Z4)/2  } 

=  {  (Z 1 1 +Z 1  )(Z  1 1  +Z2)-Z  1 22  }  {  (Z33+Z3)(Z33+Z4)-Z342  ) 

-  4  Z132  {  Z 1 1  -Z 1 2+(Z  1  +Z2)/2  }  {  Z33-Z34+(Z3+Z4)/2  } 
where  the  following  matrix  relation  is  used  : 


all  al2  00 

a21  a22  c  0  =  CAHB)  -  c2all  b22 

0  c  bl  1  bl2 
0  0  b21b22. 


B.  EIGENVECTOR  FOR  AN  ODD  OR  EVEN  RESONANCE 


Insert  the  odd  resonance  condition 

Z1 1  -  Z12  +  ZI  =  0  (B-l) 

into  Eq.  (2-3)  to  obtain 

Z12  ( II  +  12  )  +  Z13  ( 13  +  14)  =0  (B-2a) 

Z12  ( II  +  12  )  +  Z13  ( 13  +  14)  =0  (B-2b) 

Z13  ( II  +  12  )  +  (  Z33  +  Z3  )  13  +  Z34 14  =  0  (B-2c) 

Z13  ( II  +  12 )  +  Z34 13  +  ( Z33  +  Z3  )  14  =  0  .  (B-2d) 

Subtracting  Eq.  (B-2d)  from  Eq.  (B-2c),  we  obtain 

( Z33  +  Z3  -Z34 )  13  -  ( Z33  +  Z3  -  Z34 )  14  =  0 
or  13  =  14.  (B-3) 

Substitution  of  Eq.  (B-3)  into  Eqs.  (B-2)  yields 

Z12  2Z13  1  11+12  * 

=  0  .  (B-4) 

2Z13  Z33+Z34+Z3  J  L  13 


As  Eq.  (B-4)  is  a  homogeneous  system  and  is  valid  for  the  arbitrary  value  of  Z3,  the 

current  vector  is  identical  to  zero,  i.e., 

11+12=0 
13  =  0 

hence,  1 1  =  - 1 2  (B-5a) 

I  3  =  ( 1 4  )  =  0  .  (B-5b) 

Similarly,  the  eigenvector  for  another  odd  resonance  is  obtained  as 

I  3  =  - 1 4  (B-6a) 

I I  =  ( 1 2  )  =  0  .  (B-6b) 

The  eigenvector  for  the  even  resonance  may  be  obtained  from  Eq.  (2-3)  as  Zl  =  Z2 


(  Z1 1  +  Z1  )  II  +  Z12  12  +  Z13  ( 13  +  14  )  =  0  (B-7a) 

Z12  II  +  ( Z1 1  +  Z1  )  12  +  Z13  ( 13  +  14 )  =  0  (B-7b) 

Z13  ( II  +  12  )  +  (  Z33  +  Z3  )  13  +  Z34  14  =  0  (B-7c) 

Z13  ( II  +  12  )  +  Z34 13  +  (  Z33  +  Z3  )  14  =  0  (B-7d) 

Subtracting  Eq.  (B-7b)  from  Eq.  (B-7a)  and  so  Eq.  (B-7d)  from  Eq.  (B-7c),  we  get 

(Zll +Z1 -Z12)(I1 -I2)  =  0  (B-8a) 

( Z33  +  Z3  -  Z34 )  ( 13  - 14  )  =  0  .  (B-8b) 

Since  the  even  resonance  condition  does  not  include  the  odd  resonance  conditions  at  a 
time,  the  impedance  factors  in  Eqs.  (B-8)  are  not  zero.  Thus, 

11  =  12  (B-9a) 

13  =  14  .  (B-9b) 


C.  DERIVATION  OF  EQS.  (2-17)  AND  (2-18) 


t 

VI 


|  ZP  ^  ic^ 

1  ZPial  i 


Zq  13 


t 

V2 

^77. V : 


T 


z 


^ib 


- 1*5  ^ 

Zq  T  V3 
V4 

77, 


Fig.  C-l  Equivalent  circuit  with  V  I  definitions 


Give  the  voltage  and  current  definition  references  as  in  Fig.  C-l  to  obtain  circuit 


equations  : 

II  +I2  =  ia  +  ic  (C-l  a) 

13  +  14  =  ia  -  ic  (C-lb) 

Za  ia  =  Zb  ib  +  Zc  ic  (C-lc) 

VI  =  Zp  II  +  Za  ia  (C-ld) 

V2  =  Zp  12  +  Za  ia  (C-le) 

V3  =  Zq  13  +  Zb  ib  (C-lf) 

V4  =  Zq  14  +  Zb  ib  (C-lg) 

From  Eq.  (C-la)  and  (C-lb),  we  get 


ia  +  ib  =  Il  +  I2  +  I3  +  I4  . 


From  Eq.  (C-lb)  and  (C-lc),  we  get 


Za  ia  -  ( Zb  +  Zc )  ib  =  -  Zc  ( 13  + 14 ) 


(C-2) 


(C-3) 


From  Eq.  (C-2)  and  (C-3),  ia  and  ib  are  obtaind : 


ia  =  4fr±4g—  ii  +  M 2±zg  .  12  + - £2 — 13  + - £2 — 14  (C-4a) 

Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc 


ib  =  — Za — 11  + - Za — 12  +  Zc+Za  13  +  .Zc+Za  14 

Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc 


(C-4b) 


Inserting  Eqs.  (C-4)  into  Eqs.  from  (C-ld)  to  (C-lg),  system  of  V-I  equations  are  obtained 


as  follows : 


Za(Zb+Zc)  \  1 1  +  Za(Zb+Zc)  12 ^ _ ZaZb _ 13^ _ Za  Zb _ 14  (C-Sa) 

Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc 


Za+Zb+Zc 


Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc 


14  (C-5b) 


Za+Zb+Zc  Za+Zb+Zc 


Za+Zb+Zc 


Za+Zb+Zc 


14  (C-5c) 


V4=  -IIh _ CA-LQ. _ 12+  LOi&C+ca) 

Za+Zb+Zc  Za+Zb+Zc  Za+Zb+Zc 


Zb(Zc+Za))  14  (C-5d) 

Za+Zb+Zc 


ZII  Z12Z13 Z13 


Z12Z11  Z13Z13 
Z13Z13Z33Z34 


(C-5e) 


Z13Z13Z34  Z33 


where 


Zll  =Zp  + 


Za  +  Zb  +  Zc 


Z12  = 


Za  +  Zb  +  Zc 


Z13  = 


Za  +  Zb  +  Zc 


Z33  =  Zq  + 


Za  +  Zb  +  Zc 


(C-6a) 

(C-6b) 

(C-6c) 

(C-6d) 


s 


Za  +  Zb  +  Zc 

Next,  we  derive  Za,  Zb,--,  Zq  as  functons  of  Z  parameters.  For  this,  we  solve 
Eqs.  (C-6).  Divide  Eq.  (C-6b)  by  Eq.  (C-6c)  to  obtain 

Z12  =  l  +  Zc  (C-7) 

Z13  Zb  . 

Divide  Eq.  (C-6e)  by  Eq.  (C-6c)  to  obtain 

Z24  =  i  +  Zc  (C-8) 

Z13  Za  .  v  } 


Modify  Eq.  (C-6c)  to  get 

1  +-Zc  =  za(— 1 - 1- 

Zb  Z13  Zb 

1  +  .ZC  =  Zb  ( — 1 - L- 

Za  Z13  Za 

Substitute  Eqs.  (C-9)  into  Eq.  (C-7)  and  (C-8)  to  obtain 


Z12  =  za(— -  -i-) 
Z13  Z13  Zb 


ZM  -  zb  ( — 1 - L. )  or  -L.  =  Zll  ( _I - L_ 

Z13  Z13  Za  Zb  Z34  Z13  Za 

Eliminate  Zb  in  Eq.  (C-lOa)  by  substituting  Eq.  (C-lOb)  to  obtain 

_  Z132  -  Z12  Z34 

Za  - - 

Z13  -  Z34 

Inserting  the  result  of  Eq.  (C- 11)  backward,  all  parameters  arc  evaluated  : 

„  Z132  -  Z12  Z34 


Zc  =  - 


Z13-Z12 
Z132  -  Z12  Z34 


Zp  =  Zll  -Z12 
Zq  =  Z33  -  Z34 


D.  NORMALIZED  IMPEDANCE  MATRIX 

For  an  N-port  network,  the  impedance  matrix  is  defined  as 


(C-9a) 

(C-9b) 


(C-lOa) 

(C-IOb) 


(C-ll) 


(C-12) 


(C-13) 


(C-14) 

(C-15) 


where 


[  V  ]  =  VI 
V2 


[I]=  II 
12 


Zll  Z12  -  Z1N 
Z21  Z22 

[Z]=  :  \  : 

ZN1  -  -  ZNN 

Now,  we  introduce  a  normalized  voltage  vector  and  a  normalized  current  vector  defined  i 

[  v  ]  =diag  [  VZof1  ]  [  V  ]  (D-2a) 

[  i  ]  =  diag  [  VZoi  ]  [  I  ]  (D-2b) 

where  Zo i  is  the  characteristic  impdance  at  the  i'-th  port. 

Substitute  Eqs.  (D-2)  into  Eq.  (D-l)  to  obtain 

[  v ]  =  diag  [ VZo? J'1  ( Z ]  diag  [  VZo7 [i] 


where 


-UHM 

[  z  ] :  normalized  impedance  matrix  with  each  element  as 


P-3) 


P-4) 


*  diag  [  a," 1  ]  =  diag  [  a,  ]*1  P-5) 

deanwhile,  voltages  and  currents  are  decomposed  into  forward  wave  component  and 
lackward  wave  component : 

[  v  ]  =  [  v  ]+  +  [  v  ]'  P-6a) 

[i]=[i]+-[*]‘  P-6b) 

vhere  superscript  (+)  refers  to  forward  quantity  and  (')  refers  to  backward  quantity.  S 
natrix  is  defined  commonly  in  terms  of  normalized  voltage  vectors  as 

[  v  ]'  =  [  S  ]  [  v  ]+  .  P-7) 

lubstitute  Eqs.  P-6)  into  Eq.  (D-3)  and  make  use  of  the  relations  that  [v]+=  [i]+  and 
v]*=  [ij~  to  obtain 


[v]+  =  [[z]  +  [U3r1[[z]-[U]][v]- 
[S]  =  [[z]  +  [U]]-1[[z]-[U]] 


P-8) 


where  [  U  ]  is  a  unit  matrix.  S  matrix,  i.e.,  S  parameters  are  calculated  regarding  to 
normalized  z  parameters  as  follows : 

[[z]  +  [U]-1[[z]-[U]]=  adiiifyl+JJJlii  [[Z]-[U]] 

det  II  [  z  ]  +  [  U  ]  II 


i  all  al2  al3  al3  zll-1  zl2  zl3  zl3 

-  al2  all  al3  al3  zl2  zll-1  zl3  zl3 

detll[z]+[U]ll  al3  al3  a33  a34  zl3  zl3  z33-l  z34 

al3  al3  a34  a33  zl3  zl3  z34  z33-l 


detll[z]+[U]ll 


ell  cl2  cl3  cl3 
cl2  ell  cl3  cl3 
cl3  cl3  c33  c34 
cl3cl3c34  c33 


where 


all  =  (z33-z34+l)  { (zll+1)  (z33+z34+l)  -  2  zl32  ) 

al2  =  -  (z33-z34+l)  {  zl2  (z33+z34+l)  -  2  zl32  } 

al3  «  -  zl3  (zll-zl2+l)  (z33-z34+l) 

a33=  (zll-zl2+l)  {  (z33+l)  (zl  l+zl2+l)  -  2  zl32  } 

a34  =  -  (zl  1-Z12+1)  (  z34  (zl  l+zl2+l)  -  2  zl32  } 

cl  1  =  al  1  (zll-1)  +  al2  zl2  +  2  al3  zl3 

cl2  =  all  zl2  +  al2  (zll-1)  +  2  al3  zl3 

cl3  =  (al l+al2)  zl3  +  al3  (z33+z34-l) 

c33  =  2  al3  zl3  +  a33  (z33-l)  +  a34  z34 

c34  =  2  al3  zl3  +  a33  z34  +  a34  (z33-l) 

det  II  [  z  ]  +  [  U  ]  II 

=  (zll-zl2+l)(z33-z34+l)  { (zll+zl2+l)(z33+z34+l)-4  zl32  } 


P-lOa) 


p-lOb) 

P-lOc) 

P-lOd) 

P-lOe) 

P-lOf) 

P-lOg) 

P-lOh) 

(D-lOi) 

P-10j) 

P-lOk) 

(D-101) 


E.  DERIVATION  OF  EQS.(3- 10) 

For  the  notational  convenience,  the  subscripts  mn  are  mostly  omitted  in  this  section 
and  some  abbreviation  forms  are  listed  in  advance : 


/«■/ 


B«  =  8mn«  i  =  1, 2, 3 

A®  =  Amn®  *  =  1,2 

B«  =  Bmn«  /=  1,2 

D,  E,  F,  G  =  Dmn,  Emn,  Fmn,  Gmn 
V,(P  =  Vmn,<Pmn 

=  ^mn.v^  i  =  1,2 

®v(l)  =  ®mn,v(^  /  =  1,2 

Sj,  S2,  S3  =  sin  B(1)hl,  sin  B(2)h2,  sin  B(3H 
Cj,  C2,  C3  =  cos  B(1)hl,  cos  B(2)h2,  cos  B(3)t 
Ctj,  Ct2»  Q3  =  Cj/Sj,  C2/S2,  C3/S3 
zA=jcoii0 

Apply  the  boundary  condition  of  Eq.  (3-7a),  i.e., 

Etl  =  Et3 

to  Eq.  (3-4a)  and  Eq.  (3-6a)  at  z  =  0,  and  take  inner  products  over  S  with  the  orthogonal 
functions  Vt  V  x  z  and  Vt  respectively  to  obtain 

f  Vt  V  x  z  •  [  Etj  -  Et3  ]  ds  =  0  at  z  =  0 

s 

/  Vt  •  [  Etj  -  Et3  ]  ds  =  0  at  z  =  0 

s 

-A<1>S1=DS3  (E-l) 

-B(1>(B(1V1)S1=G(B(3VyA3)S3  .  (E-2) 

Take  the  inner  products  over  S  with  Vt  V,  Vt  ^  x  z  for  the  boundary  condition  of  Eq. 
(3-7c)  respectively  to  obtain 


I  Vt  V .  [  Ht3  -  Htj  ]  ds  =  J  Vt  V  Pv(1)  *v(,)  ds 
s  si 

/  Vt  x  z  •  [  Ht3  -  Htj  ]  ds  =  S  Vt<Pxz*I  Pv(1)  Xv(1) 
s  si 


at  z  =  0 


at  z  =  0 


where 


{  C+  D  C3  }  (P^/z*)  -  A<*>  (P<1VzA)  C{  =  I  Pv(‘>  •  Sv*1) 
F  +  G  C3  -  B<J)  Cj  =  S  Pv*1*  •  ev<l> 


(E-3) 


5v(1)  =  J  Xv^^VtVds 

si 

=  /  Xv*1)  •  Vt<P  x  z  ds 
si 

Note  that  the  integral  region  for  ^v(1),  ®v(1)  is  reduced  to  SI. 

In  order  to  apply  the  boundary  condition  of  Eq.  (3-7b),  we  may  take  the  inner 
product  with  function  tyi  so  as  to  expect  the  same  expressions  as  ^|i,  for  convenience. 
The  H  field  function  Xp.,  however,  is  orthogonal  to  Et  field.  Hence,  we  use  Xp  x  z 
instead  which  is  parallel  to  Et  and  integration  is  done  only  over  SI  : 


|Xp  x  z  •  Etj  ds  =  0 
si 


at  z  =  0 


E  [  A(1) !  Xp(1>  x  z  »VtV  x  z  ds  +  B(1)(P(1)/yAj)  J Xp(1)  x  z ‘Vt^  ds  ]  Sj  =  0 
mn  si  si 

MN 

or  I  [  A<*>  Sp(1)  -  B(I)(P(1)/yA1)  8p(1)  ]  Sj  =  0  ,  p  =  1,  2,  3,  -  (E-5) 

mn 

where  following  vector  identities  are  used  : 

X  X  Z  *VtV  X  Z  as  X  »VtV 

X  xz’Vt^-X.Vt^xz 

Similar  process  is  taken  for  the  boundary  conditions  at  the  interface  (II)  to  extract 
unknown  coefficents  in  Es.  (3-5)  and  (3-6). 

Take  inner  products  over  S  with  Vt  V  x  z  and  Vt  respectively  for  the  boundary 
condition  of  Eq.  (3-8a)  to  obtain 


I  Vt  v  X  z  •  [  Et2  -  Et3  ]  ds  =  0 
s 

/  Vt  9  •  {  Et2  -  Et3  ]  ds  =  0 


at  z  =  t 


at  z  =  t 


**.•**«,"  »  .'ll  .*«  ,>  '»  fcS  .  "V  _■*»  *  *-  *.  .%  "•  A  * ,  \ 


i 


A<2>  S2  =  -  c  s3  (E-6) 

B<2>  ( B<2VyA2 )  S2  =  -  F  ( B(3)/yA3 )  S3  .  (E-7) 

Take  the  inner  products  over  S  with  Vt  V,  Vt  x  z  for  the  boundary  condition  of  Eq. 
(3-8c)  respectively  to  obtain 

/  Vt  V •  [  Ht2  -  Ht3  ]  ds  =  J  Vt  V Pv(2) Xv(2) ds  atz  =  t 

s  s2 

/  Vt  x  z  •  [ Ht2 - Ht3  ] ds  =  /  Vt  9  x  z  •  £ Pv(2) Xv<2)  at z  =  t 

s  s2 

A<2>  (p<2>/zA)  C2  -  {  C  C3+  D  }  (P<1VzA)  =  I  Pv(2)  •  V2>  (E-8) 

B<2>  C2  -  {  F  C3+  G  }  =  I  Pv®  •  0v®  (E-9) 

where 

4v(2)  =  /  Xv(2)  •  VtV  ds 
s2 

0v(2)  =  /  Xv{2)  •  Vt*P  x  z  ds 
s2 

Take  the  inner  pnxiuct  over  S2  with  X^ixz  for  the  boundary  condition  of  Eq.  (3-8b)  to 
obtain 

/  X^  x  z  •  Et2  ds  =  0  at  z  =  t 

s2 

MN 

£  [  A(2)  /  Xp.(2)  x  z  *VtV  x  z  ds  +  B<2)(p(2)/yA2)  J  x  z  •Vt*P  ds  ]  S2  =  0 
mn  s2  s2 


MN 

or  I  [  A<2>  V(2)  -  B(2)(P(2)/yA2>  V2)  ]  S2  =  0  ,  p  =  1 ,  2,  3,  -  (E- 10) 

mn 

Eliminate  Cmn,  Dmn,  Fmn,  and  Gmn  in  Eqs.  (E-3),  (E-4),  (E-8),  and  (E-9)  by 
substitution  of  Eqs.  (E-l),  (E-2),  (E-6),  and  (E-7)  to  obtain 

-[p(1>C1S3+p(3>S1C3]A(1)-p®S2A®  =  S3za  I  Pv^vW  (E*1  la> 

p(3)S1A(1)+[p(2)C2S3+P(3)S2C3]A(2)  =S3za  I  Pv®^v®  (E-l  lb) 

-[P(3)C1S3+(e3/Ei)p(3)S1C3]B(1)-(e3/E2)P(2)S2B(2)  =  p(3)S3I  PvO).0vO)  (E-l  lc) 
(E3/E1)P(,)S1B(1>+[p<3>C2S3+(£3/E2)p(2>S2C3IB®  =  P(3)S3I  Pv<2)»9v(2) .  (E-l  Id) 
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Solve  Eqs.  (E- 11)  for  Amn(1),  Amn(2\  Bmn(1),  Bmn(2)  to  obtain 


A(1)  =  a(H)  X  +  a(12)  X  Pv^*^v^  (E-32a) 

A(2)  =  a(2i)  £  pvtD^vd)  +  a(22)  I  Pv(2)4v(2)  (E-12b) 

Bd)  =  bdD  X  Pva).0va)  +  b02)  I  Pv<2>«ev<2>  (E-12c) 

B(2)  =  b(2i)  x  Pv(n.M>  +  b<22>  I  Pv<2>*ev<2>  (E-12d) 

where 

amn(11)  =  Ka  (  1/Sj  )  (  8<2>Ct2  +  B<3>Ct3  )  (3-12e) 

amn<12)=  Ka(  1/(S1S3))B(3)  (3-12f) 

amn<21)  =  -  Ka  ( 1/  (S2S3) )  B<3>  (3-12g) 

amn(22>  =  -  Ka  ( 1/S2  )  (  B*1*^  +  B<3>Ct3  )  (3-12h) 

Ka  =  z*  /  { 8<3>2-B<3>Ct3(B(2>Ct2+8< 1  >Ct  1)-B<1  W2>Ct  Y  Ct2  } 
bmn<n>  =  Kb(  1/Sj  )  (B<3)Ct2+(e3/e2)B(2>Ct3)  (3-12i) 

bmn<12>  =  Kb  (  1/(S1S3)  )  (£3/62)  8(2)  (3-12j) 

bmn<21>  =  -  Kb  (  1/  (S2S3) )  (£3/8!)  B(1)  (3-12k) 

bmn^  =  -  Kb  (  1/S2  )  (B<3>Ct  j-Keyei^1^)  (3-121) 


Kb  =  B<3>/  [  (eye i  )(eye2)B(1)B(2)-B(3)Ct3  { (£3/6!  )B(1)Q2+(e3/e2)B(2)Ct { }  -B^CtjC^  ] 


Finally,  substitute  Eqs.  (E-12)  into  Eq.  (E-5)  and  Eq.  (E-10)  to  obtain 

ViM  v^M 

X  Pv(1)  U|iv(11)  +  X  Pv(2)  U*iv<12>  =  0 

M  ^2  n=  1.2, 3, 

ViM  ViM 

X  Pv(1  Ujiv(21)  +  X  Pv<2>  Unv<22>  =  0 
“l  ~2  |i  =  1, 2,  3,  v2M 


(E  13a) 


(E- 1 3b) 


where 


Unv('»=  X  S,[  &{lj)  V°  ^  -  (P(0/yA,)  b(‘>3  V°  ev01  ]  (E- 1 3c) 


i,  j  =  1,  2 


. '  -  -  ,  r.  /,  V,  /  4TV  >  .  V  Wv'-/'  ■' 


V  Si 


r 
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Eqs.  (E-13)  are  refered  to  as  a  homogeneous  system  of  equations  in  terms  of  unknown 


coefficients  Pv(i)  ( i  =  1,2  ),  and  are  expressed  in  the  matrix  form  as  follows  : 

[U][P]=0 

or 


v  =  1  2 -  vlM  v=  1  —  v2M 


(E-14) 


Integral  values  of  ^vw,  are  derived  as  follows. 


First,  we  get  the  scalar  potantial 


vector  functions  as 


where 


VtV=  -  P  (  x  212  sin  212x  cos  H2y  +  y  22  COs  222x  sin  22y  ) 
a  a  b  b  a  b 


x  112  sin 


P  =  -J  Sm  5n  1 

a  b  kc 


8/  = 


1  i  =  0 

2  i*0 


(E-15a) 
(E- 15b) 


kc2  =  ( -H12-  )2  +  (  -ML  )2 
a  b 


Take  integrals  of  the  product  of  Eqs.  (E-15)  and  Eqs.  (3-13),  (3-14),  we  get 


v.vvv.. 


v.v.  - 


^x,mn,rs(1)  =  J  Xx,rs(1)»  VtVmn  ds 
si 


where 


a  (b±wl) 

=  -P  -ma-  J  sin  -S^x  sin  -XS-x  dx  /  2  cos  -Q2-y  sin 

a  0  a  a  (b-wl\  b  2 

2 

=  -  P  JJ2L  Srm  HI(s,n,wl,b)  (E 

6ij  =  {  0  i*j 

1  »=j 


HI(s,n,wl,b) 


(sJt/wl)z  -(nji/b) 


b  wl 


'■5y,mn,rs(1)  =  1  Xy,rs(1)«  Vt^mn  ds 
si 


a  HajSdj  cos  _Sg_{y-ft>-wl)  } 

=  -P  -QZL  /  cos  x  cos  -XZLx  dx  J  sin  -HU-y - - 2 _ 

b  0  a  a  X)  bJ  y ! .  { y-b/2  \2 

2  1  wl/2  ’ 

=  -  P  5rm  SI(s,n,wl,b)  (E 

where 

SI(s,n,wl,b) 

=  0  n  +  s  :  even 

{  (lbS-1..)  (n+s-1) 

=  ^f1  f  (-D2Jo(l-GJ-  -  1)  +  (-l)2Jo(l-SJ-f  |)  ]  n  +  s  :  odd  (E 

^  u  wl  L  b  wl  2 

Jo  :  Zero  order  Bessel  function  of  the  first  kind 


*x,mn.rs(1)  =  /  Xx,rs(1)»  Vt^mnxzds 
si 


(b+wl) 

r  2 


=  P  4?"  J  sin  -ffi^-x  sin  -^x  dx  /  cos  -JllLy  sin  -SS-fy  b~w^)  dy 
b  0  a  3  rb.wn  b 3  wF  2  '  y 

2 

=  P-°?-  4-  5rm  HT(s,n,wl,b)  (E 


Jyfmn  js(1)  =  J  Xy,rs(1)»  Vt^mn  x  z  ds 
si 

a  Cb+w  cos  -S%.(y  (b~wl)  ) 

=  -P  -m^  j  cos  x  cos  -21-x  dx  J  sin  -EJ-y - — - 2 -  dy 

a  0  a  a  (b-wD  b  yll-i^hO.} 2 

9  wl/2 


_  .  p  _mZL  -JL_  Smi  SI(s,n,wl,b) 
a  8m 


(E-21) 


5x,mn,rs(2)  =  J  Xx,rs(2)*  VtVmn  ds 
s2 

b  (a+vv2) 

=  -  P  /  cos  -E^-y  cos  -^-y  dy  J  ^  sin  -E^x - - - - - dx 

a  0  b  b  (a-w2)  a  Vl-{^g-}2 

= .  p  JD2L  — b_  8sn  SI(r,m,w2,a)  (E-22) 

a  on 


cos  -IS-{x  } 


=y,mn,rs(2)  =  J  5Cy,rs(2)*  VtVmn  ds 
s2 

b  (a+w2) 

=  -  P  J  sin  -Q2-y  sin  ^-y  dy  J  ^  cos  -S^x  sjn  -TZL/X-  a-.\y 2>  dx 
b  0  b  b  fa-w2)  a  2 

.  2 

=  .  P  JUL  i  5sn  HI(r,m,w2,a)  (E-23) 

b  2 


*x,mn,rs(2)  =  J  5Cx,rs(2)»  Vt^mn  x  z  ds 
s2 


(a+w21 

2 


COS  -^-{x-^^--2^  } 


-  P  -&JI-  J  cos  -^J-y  cos  -SJ-y  dy  /  sin  -®ZLx - ^2_ — .2 


0 


=  P  — -J2-  5sn  SI(r,m,w2,a) 
b  on 


(a-w2) 
2 


a  V  1  -  (  x-a/2  I2 

1  w2/2  ' 


dx 


(E-24) 


'y.mnjs^5  =  I  $^>rs*2)*  Vt^mn  x  z  ds 

b  (a+w2) 

=  -  P  -£a^-  /  sin  -^y  sin  -^-y  dy  1  ^  cos  -^^-x  sin  -HL/X  3-w2)  dx 
a  0  b  b  (a-w21  a  W2  2  7 

_  .  p  JH2L  _tL  §sn  HI(r,m,w2,a) 

E  4 


(E-25) 


Here  are  summary  of  ®v(i) : 


^x(1)  =  -  Pmn  8rm  HI(s,n,wl,b) 

a  2 

=  -  Pmn  -22 — SL.  5rm  SI(s,n,wl,b) 
b  8m 

=  -  Pmn  -212 — k-  8sn  SI(r,m,w2,a) 
a  8n 

^y(2)  =  -  Pmn  -22.  -k-  8sn  HI(r,m,w2,a) 
b  2 

®x(1)=  pmn-22-  JL  8rm HI(s,n,wl,b) 
b  2 

9y0)  = .  Pmn  -212 — a_  8rm  SI(s,n,wl,b) 
a  8m  . 

®x(2)  =  Pmn  -22 — k-  gsn  SI(r,m,w2,a) 
b  8n 

®y(2)  =  -  Pmn  -212.  Jl  5sn  HI(r,m,w2,a) 
a  2 

8j  =  f  1  i  =  0 

2  i  #  0 


HI(s,n,wl,b) 


: _ tsrown  rCosn7ttp~wtl  (  irscocnmo+w11  f< 

(S7t/wl)2  -(n7t/b)2  2b  2b 


SI(s,n,wl,b) 

{“  °  (n-SiD  In±5=Jl 

=  BXl  [  (-1)2Jo(M32-  .  -S2L)M  I)  +  (-l)2Jo(|.n2.4..S2-)ffil  I)  ] 
4  b  wl  2  b  wl  2 


V'/' 


F.  CHOICE  OF  rs  COMBINATION 

The  numbers  for  r  and  s  and  the  choice  of  x  or  y  component  for  and  are 
assigned  according  to  the  number  v  in  Eq.  (E-14).  An  example  of  rs  combination  is 
presented  here  ;  actually,  the  results  in  chapter  4  are  calculated  by  this  combination. 

Table  F-l  shows  the  example.  In  Table  F-l,  as  M-M  wall  case  with  v  =  1  (odd 
number),  for  instance,  X  (1>  is  chosen  at  interface  (I)  for  and  calculations  so  that 
the  first  combination  of  r'  =  0  and  s'  =  0  is  assigned.  Next  r'  s'  combination  for  v  =  3  is 
that  r'  =  1  and  s'  =  0.  With  v  for  even  number,  Xx  and  is  chosen. 


Table  F-l  rs-table 


Wall 

(A) 

(B) 

M 

M 

E 

M 

(1 

E 

at  Interfaced) 

v*even 

Xx  (Jy) 

v=odd 
xy  (Jx) 

r=0,1,2,~ 

S  -  1 

r'=0, 1 ,2, — 
s'=0 

r=2r'+ 1 
s=2s' 

r  =  1,2,3,-- 
s'=  1 

r'=0, 1 ,2, —  * 
s'=0 

r=2r' 

s=2s' 

r'=0,1,2,- 

s'=0 

r*=0, 1 ,2, 
s'=0 

r=2r'+ 1 
s=2s'+ 1 

Vi* 


Wall 

v=odd 

v=even 

(A) 

(B) 

Xx  (Jy) 

xy  (Jx) 

M 

M 

r'=0 

s'=0, 1 ,2, 

r'=  1 

s’=0, 1 ,2, — 

E 

M 

r‘=0 

s=0 ,1,2,-- 

r'=0 

s=0,1,2,- 

M 

E 

r‘=0 

s=0 ,1,2,-* 

r‘=  1 

s‘=  1 ,2,3, — 

r=2r' 
s=2s'+ 1 


r=2r‘+1 
s=2s  + 1 


r=2r‘ 

s=2s' 


*  The  foundamental  resonance  mode  is  described  as  r'=l(interface(I))  or 
s'=l(interface(II))  instead  of  the  first  combination  with  r'=0  or  s'=0 


Note  that  transverse  distribution  of  each  current,  i.e.,  longitudinal  or  transverse 
current,  is  represented  only  with  one  term  in  Table  F-l;  transverse  distribution  may  be  well 
expressed  only  with  one  term  while  many  terms  are  needed  for  the  expansion  of 
londitudinal  distribution  because  of  discontinuiteis. 


>  v  v  ’/v  v  v  v  v,v; •'  •.  . 


\  A  •*.  ’ 

j.  'J.'-  V  ! 
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PROGRAM  LIST 


PROGRAM  CROSS2 

★ 


* 

★ 

★ 

* 


★ 


Subroutines  in  IMSL  are  required  in  LINKing 


Real  Ea (180) 

Common  /size/a, b,  wl,  w2,  hi,  h2,  t 
Common  /Para/Freq, Erl, Er2, Er3 

Common  /Tran/MEwall,Nreso, NHmat, Nmax, icount, coef 1 
Common  /Mat  /Ea 


External  RSNTR1 

* 

OPEN (FILE- ' CROSS2 . DAT 1 , UNIT-7 , STATUS- 1  OLD ’ ) 

* 


* 


* 


write (*, 901) 


901  format  ( 


/' 

/' 

/■ 

/* 

/' 

/' 

/' 

/' 

/' 

/' 

/’ 

/* 

/' 

/' 

/' 

/* 


b/2 


Cross-Striplines 


w2 


1 

1 

1 

1- 

1 

1- 

1 

1 

1 

0- 

0 


i 

i 

i 

- i- 


i 

i 

i 

-i - 


i  x 
-i — i- 
i  i 
i  i 
i  i 


1 

1 

1 

-1 

lwl 

-1 

1 

1 

1 


a/2 


1986-0227, 1986-0404' 
by  Tom  Uwano ' / 


- h2» 

1  Er2  (II)  1 

1  w2  1  ' 

(III)  1 - *** - 1  t  ' 

Er3  1— — — — — — — — 1  0  ' 

1  1 
1  Erl  (I)  1  ' 

- hl . 


Frequency  [MHz]  '/) 


write (*,*) 'Input  0 (MMwall) , 1 (EMwall, x-a/2)  or  2 (MEwall, y-b/2) • 
read(*,*)  MEwall 
if (MEwall. EQ.0) 

-  write <*,*)' Input  0 (combination) , 1 (region (I) )  or  2 (region (II) ) ' 
if (MEwall .EQ. 1)  write (*,*)' Input  0 (combination)  or  1 (region (I) ) ' 
if (MEwall .EQ. 2)  write (*,*)' Input  0 (combination)  or  2 (region (II) ) ' 
read(*,*)  Nreso 

if ( (MEwall. EQ.0. AND. Nreso.NE. 2) . OR. MEwall. EQ. 1)  then 
write (*,*)' Input  approximate  a  in  [mm] ' 
read(*, *) a 

write (*,*)' Input  b (fixed)  in  [mm]' 
read(*,*)  b 
else 

write {*,*)' Input  approximate  b  in  [mm]’ 


read(*,*)  b 

write (*,*)' Input  a (fixed)  in  [mm] ' 
read(*, *) a 
endif 

write (*,*)' Input  hl,t,h2  in  [mm]' 
read(*, *)  hl,t,h2 
write (*,*)' Input  wl,w2  in  [mm]' 
read(*,*)  wl,w2 

write (*,*)' Input  frequency  in  [MHz] ' 
read(*,*)  freq 

write (*,*) 'Input  Erl,Er2,Er3' 
read(*,*)  Erl,Er2,Er3 

a  -  a*l.e-3 
b  -  b*l .e-3 
hl«hl*l.e-3 
t  —  t*l.e-3 
h2«h2*l.e-3 
wl«wl*l .e-3 
w2-w2*l .e-3 
f req-freq*l .e-3 

if ( (MEwall . EQ . 0 . AND .Nreso .NE . 2) .OR.MEwall.EQ. 1)  then 
ax-a 
else 
ax-b 
endif 

write (*,*)  'Input  m,n  max' 
read ( * , * )  Nmax 

write (*,*) 'Matrix  size  (2N  x  2N)  .  Input  N’ 

read(*, *)NHmat 

Np-2*NHmat-l 

q-f req*l .e3 
ME-MEwall 

write  (7,  902)  a, b,  hi,  t,  h2,  wl,  w2 , q, Erl,  Er2,  Er3, Nmax, NHmat, ME, Nreso 
902  format  ( 'a  -',E13.6/'b  -' , E13 . 6/ 'hi-' , E13 . 6 

/'t  -',E13.6/'h2-',E13.6/’wl-',E13.6/’w2-',E13.6 
/ 'freq-',E13.6,  '  MHz ' / 'Erl, Er2, Er3-' , 3F7 .2 
/ 'm, n  max  - ' , 15 

/'Matrix  2N  x  2N,  N- ', 13/ 'MEwall- 13/ 'Nreso  -’,13) 
icount-0 
coef 1-1 . 

Call  ZER02 (RSNTR1, ax, ier) 

if ( (MEwall. EQ.O. AND. Nreso. NE. 2) . OR. MEwall .EQ. 1)  then 
write(*,*)'  final  a  in  [mm]  -' , ax*l .e3, ier 
write(7,*)'  final  a  in  [mm]  - ' , ax*l .e3, ier 


else 

write(*,*)'  final  b  in  [mm]  , ax*l .e3, ier 
write(7,*)'  final  b  in  [mm]  - ' , ax*l .e3, ier 
endif 


write (*, ' (4E16.7) ')  1.,  (Ea (i) , i-l,Np) 
write (7, * (4E16.7)  •)  1.,  (Ea (i) , i-1, Np) 


FUNCTION  RSNTRl(ax) 
Cross-St riplines 


1986-0226,1986-0401 
by  Tom  Uwano 


1 

i 

i 

i 

1 

i 

i 

i 

1 

i 

i 

i 

1 - 

— i_ 

-i - 

— i 

1 

i 

i 

i 

1 - 

- i- 

-i - 

— i 

1 

i 

i 

i 

1 

i 

i 

i 

1 

i 

i 

i 

a - 

1 

1 

(III)l- 
Er3  1- 
1 
1 


«=»====!  o 
1 

(I)  1 


Frequency  [GHz] 

Integer*4  MNmax, mn 

Real  muO, KKO, Kkl, Kk2, Kk3 

Real  U(2,2) , El (32400) , Ed (180) , wk (180) , Esl (8100)  ,Es2 (8100) 
Real  Ep(32041) ,Ea(lB0) 

Complex  JJ 

Complex  zhat, yhatl, yhat2 

Complex  Btmnl , Btmn2 , Btmn3 

Complex  all, al2, a21, a22,bll,bl2,b21, b22 

Complex  SI1, SI2, SI3, COl, C02 , C03 , TA1 , TA2 , TA3 

Complex  dela, delb, zdela, b3delb 

Complex  x, FARG, CMCOS, CMS IN 

Common  /size/a, b, wl, w2, hi, h2, t 

Common  /Para/Freq, Erl,Er2,Er3 

Common  /Tran/MEwall,Nreso,NHmat,Nmax, icount,coefl 
Common  /Mat  /Ea 

Parameter  (JJ- (0 . , 1 . ) , pai-3 .1415926) 

Parameter  (epsiO-8 . 855e-12,mu0-pai*4 .e-7 ) 


Function  definition 


FARG (X) -REAL (X) +JJ*SIGN < 1 . , AIMAG (X) ) *AMIN1 (70 . , ABS (AIMAG (X) ) ) 
CMCOS (X) -CCOS (FARG (X)  ) 

CMSIN (X) — CSIN (FARG (X) ) 


if ( (MEwall .EQ. 0 . AND .Nreso .NE . 2) . OR .MEwall . EQ . 1)  then 
a-ax 


endif 

Nmat-2*NHmat 

Np«Nmat-l 

NHp-NHmat-1 

Er31=Er3/Erl 

Er32«Er3/Er2 

Omega-2 . *pai*Freq*l .e9 

KkO=Omega*Omega*muO*epsiO 

Kkl— Erl*KkO 

Kk2«Er2*KkO 

Kk3— Er3*KkO 

zhat-JJ*Oraega *mu0 

yhatl-JJ*Omega*epsi0*Erl 

yhat2«JJ*Omega*epsi0*Er2 


-  Matrix  Size 

Mumax-NHmat 

Numax-Mumax 


N  x  N  (N—  2  x  Mumax) 


Nmax-Nmax 

MNmax-Nmax*Nmax 

do  201  mu-1, Mumax 
do  201  nu-l,Numax 

Call  RSTBL(irl, isl, ixyl,mu,  1) 
Call  RSTBL (ir2, is2, ixy2,mu, 2) 
Call  RSTBL ( jrl, jsl, jxyl,  nu,  1) 
Call  RSTBL (jr2, js2, jxy2,  nu,  2 ) 


- Begining  of  Summation  in  terms  of  M,N - 

U(l,l)-0. 

U(l,2)-0. 

U(2, 1) -0 . 

U(2, 2) -0 . 

do  202  mn-1, MNmax 
Call  MNTBL (m, n,mn) 

if (irl .NE.m.AND . jrl .NE.m. AND . is2 .NE.n.AND. js2 .NE.n)  go  to  1 

GGmn- (Float (m) *pai/a) **Z r (Float (n) *pai/b) **2 
Btmn 1 -CSQRT ( CMP  LX ( Kk 1 -GGmn ) ) 

Btmn2-CSQRT (CMP LX (Kk2-GGmn) ) 

Btmn3-CSQRT (CMP LX (Kk3-GGmn) ) 

C01-CMC0S (Btmnl*hl) 

C02-CMC0S (Btmn2*h2) 

C03-CMC0S (Btmn3*t  ) 

511- CMSIN (Btmnl*hl) 

512- CMSIN (Btmn2*h2) 


S 1 3-CMS IN ( Btmn3  *t  ) 

TA1-SI1/C01 

TA2-SI2/C02 

TA3-SI3/C03 

dela-Btmn3*Btmn3*TAl*TA2*TA3-Btmn3* (Btmn2*TAl+Btmnl*TA2) 

-Btmnl*Btmn2*TA3 

de lb-Btmn 1  * Btmn2 * E r 3 1 * E r 3 2 * T A1 * TA2 * TA3 

-Btmn3* (Btmnl*Er31*TAl+Btmn2*Er32*TA2) -Btmn3*Btmn3*TA3 
zdela-zhat/dela 

all-  zdela*(Btmn2*TA3+Btmn3*TA2) 
al2«  (1 . /C03) *zdela*Btmn3*TA2 
a21— (1./C03)  *zdela*Btmn3*TAl 
a22— zdela*  (Btmnl*TA3+Btmn3*TAl) 
b3delb-Btmn3/delb 

bll-  b3delb* (Btmn3*TA3+Btmn2*Er32*TA2) 
bl2-  (1./C03) *b3delb*Btmn2 *Er32 *TA2 
b21=- (1 . /C03) *b3delb*Btmnl*Er31*TAl 
b22=-b3delb* (Btmn3*TA3+Btmnl*Er31*TAl) 

Call  GSITH1 (Gmul, Tmul, ixyl,m, n, irl, isl) 

Call  GSITH2 (Gmu2, Tmu2, ixy2,m, n, ir2, is2) 

Call  GSITH1 (Gnul, Tnul, jxyl,m, n, jrl, jsl) 

Call  GSITH2 (Gnu2, Tnu2, jxy2,m, n, jr2, js2) 

Ul 1— AIMAG (  TAl*Gmul*all*Gnul-TAl*Btmnl/yhatl*Tmul*bll*Tnul) 
U12— AIMAG (  TA1 *Gmul *al2 *Gnu2 -TA1 *Btmnl /yhat 1 *Tmul *bl2  *Tnu2 ) 
U2 1-AIMAG (  TA2 *Gmu2  *a2 1 *Gnul-TA2  *Btmn2 /yhat2  *Tmu2  *b2 1 *Tnul ) 
U22-AIMAG (  TA2 *Gmu2*a22*Gnu2 -TA2 *Btmn2 /yhat 2  *Tmu2 *b22 *Tnu2 ) 

U(l, 1)-U<1, 1)+U11 
U(l,2)-U(l,2)+U12 
U ( 2 , 1 ) -U (2,1) +U2 1 
U(2,2)-U(2,2) +U22 


1  continue 
202  continue 


End  of  Summation 


do  203  i-1,2 
do  203  j-1,2 

Call  MXTBL(k, 1, i, j , nu, mu, Numax) 
kl-k+Nmat* (1-1) 

El (kl) -U{i, j) 

203  continue 

201  continue 


- process  to  avoid  overflow  in  the  determinant 

icount-icount+1 
if (icount .EQ . 1)  then 
sum-0 . 


•w* 
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do  205  i-l,Nmat 
kl-(i-l) *Nmat+i 
3 um* 3 um+ ALOG  (ABS  (El  (kl)  )  ) 
coef 1-EXP ( (10. -sum) /FLOAT (Nmat) ) 
else 
endif 


do  206  i-l,Nmat*Nmat 
El (i) -coef 1*E1 (i) 


end  of  process 


do  204  i«l, NHmat 
do  204  j=l, NHmat 
ij=i+NHmat* ( j-1) 
kll=i+Nmat* ( j-1) 
kl2=i+NHmat+Nmat* ( j+NHmat-1) 
Esl  (i j ) =E1 (kll) 

Es2  (ij)=El(kl2) 
continue 


-Finding  Pnu — 


if (Nreso.EQ.0)  then 
do  207  i«l,Np 
Ea(i)—El(i) 
do  208  i=l,Np*Np 
kl=Nmat+i+(i-l) /Np 
Ep(i)-El(kl) 

Call  LINV3F  (Ep,  Ea,  2, Np, Np, dl, d2,  wk,  ier) 
endif 

i f (Nreso . EQ . 1 . AND . NHmat . NE . 1 )  Ea(l)»-Esl (1) /Esl (NHmat+1) 
if (Nreso. EQ. 2. AND. NHmat. NE.l)  Ea (1)=-Es2 (1) /Es2 (NHmat+1) 


Call  LINV3F (Esl, Ed, 4, NHmat, NHmat, dl, d2 , wk, ier) 
detl*=dl*2 . **d2 

Call  LINV3F (Es2, Ed, 4, NHmat , NHmat, dl , d2 , wk, ier) 
det2-dl*2.**d 2 

Call  LINV3F (El, Ed, 4, Nmat, Nmat, dl, d2, wk, ier) 

dett=dl*2.**d2 

if (Nreso. EQ.l)  RSNTRl-detl 

if (Nreso. EQ. 2)  RSNTRl-det2 

if (Nreso.EQ.0)  RSNTRl-dett 

if ( (MEwall .EQ. 0 . AND .Nreso .NE. 2) . OR.MEwall .EQ. 1)  then 
write (*, *) a*l .e3,dett,detl,det2 
write (7, *)a*l.e3,dett,detl,det2 
else 

write (*, *) b*l .e3, dett, detl, det2 
write (7, *)b*l .e3, dett, detl, det2 
endif 

return 


i r i *\mT •  4  * ~ 


v  /  v  *.vv  .V  .V 

\  V I r 


END  of  RSNTR1 


-  Begining  of  Gai  ans  Theta  calculation 

SUBROUTINE  GSITH1 (Gl, Tl, ixy, im, in, ir, is) 

Gsai  and  Theta  on  the  interface  btwn  (I)  and  (III) 

1986-0226, 1986-0402 

Real  m, n 

Common  /aize/a, b, wl, w2 
Parameter  (pai-3 . 1415926) 

DL (I) -FLOAT (MIN (1+1,2) ) 

DLIJ (I, J) “FLOAT ( 1/ (IABS (I-J)+l) ) 


if (ixy .EQ. 1)  go  to  1 
if (ixy. EQ. 2)  go  to  2 

-  x  and  interface  (I) 

1  continue 
if(ir.EQ.im)  then 

pma-m*pai/a 
pnb-n*pai/b 
GGmn-pma  *pma +pnb  *pnb 
Gmn-SQRT (GGmn) 

Pmn=SQRT (DL(im) *DL(in) / (a*b) ) /Gmn 
DD=DLIJ (ir, im) 

Fan=FIl (ia, in, wl, b) 

Gl— Pmn*pma* (a/2 . ) *DD*Fan 
Tl-  Pmn*pnb* (a/2 . ) *DD*Fan 
elae 
Gl-0 . 

Tl-0. 

endif 

return 

-  y  and  interface  (I) 

2  continue 
if(ir.EQ.im)  then 

pma-m*pai/a 
pnb«n*pai/b 
GGmn-pma  *pma+pnb*pnb 
Gmn-SQRT (GGmn) 

Pmn-SQRT (DL ( im) *DL(in) / (a*b) ) /Gmn 
Dl-DL(im) 

DD-DLIJ (im, ir) 


'ft 


Fsn-FI2  (i.3,  in, wl,b) 

Gl— Pmn*pnb*  (a/Dl)  *DD*Fsn 
Tl— Pmn*pma*  (a/Dl)  *DD*Fsn 
else 
Gl-O. 

Tl-O. 

endif 

return 

end 


SUBROUTINE  GSITH2 (G2, T2, ixy, im,  in, ir, is) 


Gsai  and  Theta  on  the  interface  btwn  (II)  and  (III) 

1986-0226,1986-0402 

Real  m, n 

Common  /size/a, b, wl, w2 
Parameter  (pai«3 . 1415926) 

DL (I) -FLOAT (MIN (1+1,2) ) 

DLIJ (I, J) —FLOAT (1/ (IABS (I-J) +1) ) 


if (ixy .EQ. 1)  go  to  1 
if (ixy. EQ. 2)  go  to  2 

k 

\ - x  and  interface  (II) 

t 

1  continue 

if(is.EQ.in)  then 
pma=m*pai/a 
pnb«n*pai/b 
GGmn-pma  *pma+pnb*pnb 
Gmn-SQRT (GGmn) 

Pmn-SQRT (DL (im) *DL(in) / (a*b) > /Gmn 
Dl-DL(in) 

DD— DLIJ (is, in) 

Fsn-Fl2 (ir, im, w2, a) 

G2=-Pmn*pma* (b/Dl) *DD*Fsn 
T2—  Pmn*pnb* (b/Dl) *DD*Fsn 
else 
G2-0. 

T2-0. 

endif 

return 

-  y  and  interface  (II) 

2  continue 

if (is. EQ. in)  then 
pma-m*pai/a 


**  -  -A 


v®.*- w 


m 


pnb-n*pai/b 
GGmn-pma  *pma  +pnb  *pnb 
Gmn-SQRT (GGmn) 

Pmn«SQRT(DL(im) *DL(in) / <a*b) ) /Gmn 
DD**DLI  J (is,  in) 

Fsn-FIl (ir, im, w2, a) 

G2«-Pmn*pnb* (b/2 . ) *DD*Fsn 
T2=-Pmn*pma* (b/2 . ) *DD*Fsn 
else 
G2**0 . 

T2=0 . 
endif 
return 
end 

end  of  Gsi  and  Theta 


FUNCTION  FI1 (is, in, w, b) 

Real  n 

pai=3. 1415926 

s=is 

n-in 

psw“s*pai/w 
pnb=n*pai/b 
if (pnb.EQ.psw)  then 

Fll»-w/2.*SlN( (n-s) *pai/2 . ) 
else 

h=  psw/ (pnb*pnb-psw*psw) 

c-COS (pnb* (b-w) /2 . ) - (-1 . ) **is*COS (pnb* (b+w) /2 . ) 
Fll=-h*c 
endif 
return 
end 

FUNCTION  FI2 (is, in, w,b) 

Real  n 

pai-3. 1415926 

s»is 

n-in 

if (MOD( (in+is) ,2) .EQ.O)  then 
FI2-0. 
else 

argl-ABS (n*w*pai/ (2 . *b) +s*pai/2 . ) 
arg2-ABS (n*w*pai/ (2 . *b) -s*pai/2 . ) 
h-BSJO(argl) * (-1 . ) **(  (in+is-1) /2) 
c-BSJO (arg2) *  (-1 . )  *M (in-is-1) /2) 

Fl2*w*pai/4 . * (h+c) 
endif 
return 


*  *  r.  ~  ' 


ena 

* 

* 

SUBROUTINE  MXTBL (k, 1, i, j , Nu, Mu, Nmax) 

k=Mu+Nmax* (i-1) 

l«Nu+Nmax* ( j-1) 

return 

end 

* 

* 

SUBROUTINE  RSTBL (ir, is, inf oxy, NNu, i) 

Common  /Tran/MEwall, Nre30, NHmat,Nmax, icount, coef 1 

if(i.EQ.l)  go  to  1 
i f ( i . EQ . 2 )  go  to  2 


1  continue 

- Interface  (I) — 

if (MOD (NNu, 2) .EQ. 0)  then 

r,s  for  the  Xx(l) 

infoxy«l 

n=NNu/2 

iir-n-1 

if (MEwall.EQ.l)  iir=iir+l 
iis-1 

if (MEwall.EQ.2)  iis-0 
if (MEwall.EQ.O)  then 
ir=2*iir+l 
i3=2*iis 
endif 

if (MEwall.EQ. 1)  then 
ir-2*iir 
is-2*iis 
endif 

if (MEwall.EQ. 2)  then 
ir=2*iir+l 


NNu  seven 


is=2*iis+l 

endif 


else 


infoxy-2 
n-(NNu+l) /2 
iir-n-1 
iis-0 

if (MEwall.EQ.O)  then 
ir-2*iir+l 
is-2*iis 
endif 

if (MEwall.EQ.l)  then 
irdum-iir 

if (irdum.EQ.O)  iir-1 
if (irdum.EQ. 1)  iir-0 
ir-2*iir 
is-2*iis 


r,s  for  the  Xy(l) 


NNu : odd 


endif 

if (MEwall.EQ.2)  then 
ir=2  *iir+l 
13=2*113+1 
endif 

endif 

return 

2  continue 

- Interface  (ii) — 

if (MOD (NNu, 2 ) . EQ . 0 )  then 

r,s  for  the  Xy(2) 

infoxy=2 

n=NNu/2 

iis=n-l 

if (MEwall . EQ . 2)  iis=iis+l 
iir=l 

if (MEwall. EQ.l)  iir=0 
if (MEwall. EQ.O)  then 
ir=2*iir 
is=2*iis+l 
endif 

if (MEwall. EQ.l)  then 
ir=2*iir+l 
is=2*iis+l 
endif 

if (MEwall .EQ. 2)  then 
ir=2*iir 
is=2*iis 
endif 

else 


NNu: even 


r,s  for  the  Xx(2) 


NNu : odd 


infoxy=l 
n= (NNu+1) /2 
iis=n-l 
iir=0 

if (MEwall . EQ . 0)  then 
ir=2*iir 
is=2*iis+l 
endif 

if (MEwall. EQ. 1)  then 
ir=2*iir+l 
is=2*iis+l 
endif 

if (MEwall. EQ. 2)  then 
isduin“iis 

if (isdum.EQ. 0)  iis-1 

if (isdum.EQ. 1)  iis=0 

ir=2*iir 

is=2*iis 

endif 

endif 

return 

end 
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SUBROUTINE  MNTBL (m, n, MN) 

Common  /Tran/MEwall,Nreso,NHmat,Nmax, icount, coef 1 

M  -  M  Wall  case 
Integer*4  MN 
n=(MN-l) /Nmax 
m=MN-Nmax*n-l 
if (MEwall .EQ.O)  then 
m=*2*m+l 
n=2*n+l 
endif 

if (MEwall .EQ. 1)  then 
m=2*m 
n»2*n+l 
endif 

if (MEwall . EQ . 2)  then 
m=2*m+l 
n-2*n 
endif 
return 
end 


FUNCTION  BSJO(x) 

Bessel  function  19 86-0226 

ax=ABS (x) 

if(ax.GT.3.)  go  to  1 

x2=x*x/9 . 

x4=x2*x2 

x6=x2*x4 

x8=x4*x4 

xl0=x2*x8 

xl2«x6*x6 

BSJ0-1.-2 .2499997d0*x2+l . 2656208d0*x4-0 . 3l63866d0*x6 
+0 .044447 9d0*x8 -0.003 944 4d0*xl0+0 .0002 Id0*xl2 

return 
1  continue 
xl=3 . /ax 
x2=xl*xl 
x3«xl*x2 
x4=xl*x3 
x5-xl*x4 
x6»xl*x5 

f-  0.797884s6d0-0.00000077d0*xl-0 . 0055274d0*x2-0 . 00009512d0*x3 
-  +0.00137237d0*x4-0.00072805d0*x5+0.00014476d0*x6 

t*ax-0 . 78539816d0-0 .04166397 d0*xl-0.00003954d0*x2 +0.0026257 3d0*x3 
-0. 00054 125d0*x4-0. 0002 9333d0*x5+0 . 00013558d0*x6 
BSJO-f *cos (t) /SQRT (ax) 
return 
end 


SUBROUTINE  ZER02 (Func, x,  ier) 


ier=l  :  no  root 

ier=2  :  dosen't  converge 

5%  step  toward  + 

Implicit  Real*8 (a-h, o-z) 
Real  x, Func, ylast 
ier=0 
xl=x 

step=0 . 05*DABS (xl) 
yl=Func (SNGL(xl)  ) 
ylast=SNGL (yl ) 

if (ylast . EQ . 0  . )  then 
x=xl 
return 
endif 

x3=xl+step 
y3=Func (SNGL(x3) ) 
ylast=SNGL (y3) 
if (ylast .EQ. 0 . )  then 
x*x3 
return 
endif 


1986-0321, 1986-0502 
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if (yl *y3 . LT . 0 . )  then 

x2=xl+ (x3-xl) *yl/  (yl-y3) 
y2=Func (SNGL(x2) ) 
xlast=x2 
yla3t=SNGL(y2) 
if (ylast .EQ. 0 . )  then 
x=xlast 
return 
endif 

else 

xnew=xl+ (x3-xl) *yl/ (yl-y3) 
del=y3-yl 

if (del*yl . LT . 0 . )  then 

forward  step 

x2»x3 

y2=y3 

x3“x3+ (x3-xl) 

if (xnew. LT . x3)  x3=xnew+0 . 01 *DABS (xnew) 

y3“Func (SNGL(x3) ) 

xlast=x3 

ylast»SNGL(y3) 

if (y2 *y3 . LE . 0 . }  go  to  100 

do  201  i-1,5 

xnew-x3+ (x3-xl ) 

xl-x2 


rJ>  'JBrawnwwwwwHwiT 
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x2-x3 

x3-xnew 

yl-y2 

y2-y3 

y3-Func(SNGL(x3) ) 

xlast«x3 

ylast«SNGL (y3) 

if (y2*y3 . LE . 0 . )  go  to  100 

201  continue 
ier«l 
x-x3 
return 
else 

*  backward  step 
x2=xl 

y2-yl 

xl=xl- (x3-xl) 

if (xnew.GT.xl)  xl-xnew-0 .01*DABS (xnew) 

yl-Func (SNGL(xl) ) 

xlast*xl 

ylast*SNGL (yl) 

if (y2*yl.LE.0.)  go  to  100 

do  202  i-1,5 

xnew-xl- (x3-xl) 

x3«x2 

x2-xl 

xl“xnew 

y3-y2 

y2*yl 

yl=Func (SNGL(xl) ) 

xlast=xl 

ylast=SNGL (yl) 

if (y2*yl . LE . 0 . )  go  to  100 

202  continue 
ier=l 
x=xl 
return 
endif 

endif 

* 

100  continue 

if (ylast .EQ. 0 . )  then 
x**xlast 
return 
endif 

* 

do  203  i-1,10 
if(i.NE.l)  then 
y2-Func(SNGL(x2) ) 
xlast«x2 
ylast-SNGL(y2) 
if (ylast .EQ. 0 . )  then 
x-xlast 
return 
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endif 

endif 

Call  ABCNEW(xnew, xl, x2, x3, yl, y2, y3) 
if (yl*y2 .GT . 0 . )  then 
xl-x2 
yl-y2 
else 
x3“x2 
y3-y2 
endif 
x2«xnew 
epsi-l .d-5 

dela-DABS ( (x2-xlast) /x2) 
if (dela . LT .epsi)  then 
x-x2 
return 
endif 

203  continue 
ier-2 
x-x2 
return 
end 


SUBROUTINE  ABCNEW (xnew, xl , x2 , x3, yl , y2 , y3 ) 
Implicit  Real*8(a-z) 

Call  ABC(a,b,c,xl,x2,x3,yl,y2,y3) 
bb-b 

if (bb . EQ . 0 . )  bb-l.d-30 
if (DABS(a/bb) .LT.l.e-10)  then 
xnew— c/b 
else 

ac-b*b-4 . *a*c 
anxl- (-b-DSQRT (ac) ) / (2 . *a) 
anx2- ( -b+DSQRT (ac) ) / (2.*a) 
xnew-anxl 

i f (anx2 . GE . xl . AND . anx2 . LE . x3)  xnew*anx2 
endif 


SUBROUTINE  ABC (a, b, c, xl , x2 , x3 , y 1 , y2 , y3 ) 

Implicit  Real*8(a-z) 

all—  1 .  /  ( (x2-xl>  *  <xl-x3) ) 

a  12—1 .  /  ( (x2-xl)  *  <x3-x2) ) 

a  13—  1 .  /  <  (x3-x2)  *  (xl-x3) ) 

a21--all» (x3+x2) 

a22— al2*  (xl+x3) 

a23— al3*  (x2+xl) 

a31-all*x2*x3 

a32-al2*xl*x3 

a33-al3*xl»x2 

a-all*yl+al2*y2+al3*y3 


I'ssasaaag  j  wKWLf 


b-a21*yl+a22*y2+a23*y3 

c-a31*yl+a32*y2+a33*y3 

return 

end 


